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(57) Abstract 

A system is provided for performing real-time warping in one or two dimensions (e.g. dynamic zoom and shrink in both horizontal 
and vertical directions), and temporal warping of image frames (e.g. fast motion and slow motion). The system accomplishes warping 
with proper Nyquist bandlimiting. Three different processes are performed in sequence to create the warped image according to the present 
invention. The first process maps a desired target pixel or frame spacing to an appropriate FIR filter wherein the size of the FIR filter is 
chosen based on the desired target pixel or frame spacing. Once the filter is chosen, a second process creates an appropriate number of 
upsamplcd pixels, lines or frames based on the input source pixels, lines or frames. The last process multiplies the correct filter coefficients 
from the first process with the interpolated pixels, lines or frames from the second process, resulting in correctly Nyquist bandlimited target 
pixels, lines or frames. 
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1 

DIGITAL IMAGE WARPING SYSTEM 

Field of the Invention 

5 The present invention relates in general to digital 

signal processing, and more particularly to an image 
warping system for generating an output stream of warped 
pixels and/or lines and/or frames from an input stream of 
pixels and/or lines and/or frames, respectively. 

10 

Background of the Invention 

Digital image warping is a process of dynamically 
resampling a regularly spaced source image to produce a 

15 target image with different spacing than the source 

image. Many modern movies use digital warping to create 
mind boggling effects such as perspectives, odd shaped 
image boundaries, and texture mapping, to name but a few. 
Prior art real time image warping systems are extremely 

20 complicated and expensive due to the massive number of 

computations which must be performed in real time. Prior 
art systems do not dynamically bandlimit the spatial 
frequency of the warped image to the Nyquist limit, 
causing unwanted aliasing distortion. Examples of well 

25 known prior art digital image warping systems are 

described in the following publications and patent: 

I. Wolberg, George, "Separable Image Warping: 
Implications and Techniques", Ph.D. Thesis, Dept. of 
Computer Science, Columbia University, NY, 1990, and 

30 "Digital Image Warping", IEEE Computer Society Press, Los 

Alamitos, CA, 1990; 

II. U.S. Patent 5,355,328 (Arbeiter, et al) . 
Summary of the Invention 

35 

According to the present invention, an image warping 
system is provided for performing real-time dynamic zoom 
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and shrink in both horizontal and vertical directions 
with proper Nyquist bandlimiting, and requiring no 
external line stores or other memory. This same 
invention can also temporally warp images (eg. dynamic 
5 changes between fast motion and slow motion) by utilizing 

external frame stores in place of registers. Accordingly, 
the system of the present invention constitutes a 
significant advance over the known prior art. 
Three different processes are performed in sequence to 

10 create the warped image according to the present 

invention. The first process maps a desired target pixel 
(or line or frame) spacing to an appropriate FIR filter 
wherein the size of the FIR filter is chosen based on the 
desired target pixel (or line or frame) spacing. Once 

15 the filter is chosen, a second process creates an 

appropriate number of upsampled pixels (or lines or 
frames) based on the input source pixels (or lines or 
frames) . The last process multiplies the correct filter 
coefficients from the first process with the interpolated 

20 pixels (or lines or frames) from the second process, 

resulting in correctly Nyquist bandlimited target pixels 
(or lines or frames) . 

Thus, in accordance with one aspect of the 
25 invention, there is provided an apparatus and method for 

digitally warping an input stream of one of either pixels 
and/or lines and/or frames to produce an output stream of 
warped pixels and/or lines and/or frames in accordance 
with a user defined target increment defining the desired 
30 spacing between output pixels and/or lines and/or frames. 

In the case of digitally warped pixels and/or lines, the 
resulting image is considered to be spatially warped. In 
the case of digitally warped frames, the resulting set of 
frames is considered to be temporally warped. 

35 

In a first embodiment of the invention, the system 
comprises a filter look-up table for mapping the target 
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increment into a specific filter size, and a pair of 
registers to pipeline three successive filter sizes 
output from the filter look-up table for input to a pair 
of comparators. A first one of the comparators 
5 determines the largest filter size from the three values 

output by the filter look-up table and both registers, 
while the second comparator determines the largest filter 
size output from the filter look-up table and the first 
register. The output of the first comparator is received 

10 by an upsample increment generator which maps the filter 

size to a filter factor and then right-shifts the target 
increment by this amount to generate an upsample 
increment. The upsample increment is received by an 
interpolator which in response generates an appropriate 

15 number of evenly spaced upsampled pixel values. The 

output from the second comparator is stored in a pipeline 
comprising two further registers connected respectively 
to left and right half kernel filter generators. The 
left and right half kernel filter generators receive and 

2 0 multiply the upsampled pixels from the interpolator by 

respective half kernel coefficients to create the warped 
output pixels. 

In a second embodiment of the invention, the system 
25 comprises a filter look-up table for mapping the target 

increment into a specific filter size, and a single 
register for holding the filter size from the filter 
look-up table for input into left and right half kernel 
filter generators. The output of the filter look-up table 

3 0 is received by an upsample increment generator which maps 

the filter size to a filter factor and then right-shifts 
the target increment by this amount to generate an 
upsample increment. The upsample increment is received by 
an interpolator which in response generates an 
3 5 appropriate number of evenly spaced upsampled pixel 

values. The left and right half filter kernel generators 
receive and multiply the upsampled pixels from the 
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interpolator by respective half kernel coefficients to 
create the warped output pixels. 

Each embodiment can be extended to digitally warp 
5 lines and/or frames by replacing the registers with line 

stores and frame stores, respectively. These larger 
memory storage units perform the identical function as 
the registers in the embodiments described in detail 
below. 

10 The warping technique of the present invention may be 
implemented utilizing the image filtering and image 
extension apparatus and methods described in co-pending 
U.S. patent application numbers 08/126,388; 08/125,530; 
08/172,065; 08/124,201, all assigned to Genesis Microchip 

15 Inc. , the contents of which are incorporated herein by 

reference. 

Furthermore, the digital image warping system 
according to the present invention builds on technology 

20 described in U.S. Patent 5,355,328 (Arbeiter et al) , 

referred to above, which relates to a method of resizing 
source images to an arbitrary but fixed target image 
size. According to this prior art system, each target 
pixel and/or line and/or frame is resized by the same 

25 amount relative to the source pixel and/or line and/or 
frame. According to the present invention, each target 
pixel and/or line and/or frame may be individually sized 
arbitrarily relative to the source pixel and/or line 
and/or frame. More specifically, according to the 

30 present invention, source pixels and/or lines and/or 

frames are dynamically upsampled based on changing target 
pixel and/or line and/or frame resize factors. The 
source pixels and/or lines and/or frames are always 
upsampled to the nearest first octave and then decimated 

35 by fixed sized FIR filters to create Nyquist bandlimited 

target pixels and/or lines and/or frames. The filter 
sizes are dynamically determined for each target pixel 
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and/or line and/or frame based on spatial positioning or 
temporal positioning as dictated by the user. Major 
silicon cost reduction is achieved with the present 
invention over prior art video warping systems through 
5 the use of fixed size FIR filters. 



prief Description of the Drawings 

A detailed description of the preferred embodiment is 
10 provided herein below with reference to the following 
drawings, in which: 

figure 1 shows a plurality of target pixels derived 
irom a plurality of source pixels , wherein the target 
15 pixels are spaced according to the user defined target 

increment; 

Figure 2 is a block diagram showing functional 
components of the digital image warping system according 
20 to the present invention; 

Figure 3 shows the relative spatial positioning of 
source and target pixels for dynamically varying user 
defined target increments , and the selection of 
25 appropriate FIR filters for generating the target pixels 

using the system of Figure 2, according to a 
representative example; 

Figure 4 is a block diagram showing functional 
3 0 components of the digital image warping system according 

to the preferred embodiment of the present invention; and 
Figure 5 shows the relative spatial positioning of the 
source and target pixels for dynamically varying user 
defined target increments, and the selection of 
35 appropriate FIR filters for generating the target pixels 
using the preferred embodiment of Figure 4. 
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Detailed Description of the Alternative and Preferred 
Embodiments 

In the description below, the preferred and 
5 alternative embodiments apply equally to the warping of 

pixels and/or lines. The same inventive principles can be 
used to temporally warp frames by substituting frames 
stores for each register in Figures 2 and 4. 
Furthermore, the invention is not limited to video image 
10 warping. For example, in the case of one-dimensional 

warping, the principles of the present invention may be 
applied to audio rather than video samples (eg. to 
correct jitter on a digital audio tape drive). 

15 As discussed above, digital warping of an image is 

specified by target pixel spacing with respect to source 
pixel spacing. The horizontal space between each target 
pixel and the vertical space between each target line 
must be specified by the user. Likewise, the temporal 

20 spacing (ie. time) between successive frames of an image 

must be user specified. This detailed information 
provides a very fine level of control that can be used to 
create very complex warps. The generation of user 
defined target spacing is known in the art, and does not 

25 form part of the present invention. The system of the 

present invention receives as inputs a stream of source 
pixels (and/or lines and/or frames) and the desired 
target increment (Tarlnc) which defines the user- 
specified spacing. 

30 

Figure 1 illustrates the generation of target pixels 
from source pixels based on a user defined target 
increment (Tarlnc) . For simplicity, the remainder of 
this disclosure will refer to the generation of target 
35 pixels from source pixels, it being understood that 

exactly the same methodology and circuitry applies to 
vertical processing by substituting video lines for 
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pixels and temporal processing by substituting video 
frames for pixels. 

As can be seen from Figure 1, the first and last 
5 target pixels of a warped image are spatially aligned 

with the first and last source pixels. However, the 
spacing of the intermediate target pixels is defined by 
the target increment (Tarlnc) . Specifically, in the 
example shown in Figure 1, the first (i.e. left-most) 
10 intermediate target pixel is spaced from the first target 

pixel by an amount equal to 1.5 times the spacing between 
the source pixels, the second intermediate target pixel 
is spaced from the first intermediate target pixel by an 
amount equal to 0.75 times the spacing between the source 
15 pixels, while the last target pixel is separated from the 

second intermediate target pixel by an amount equal to 
0.75 times the spacing between the source pixels, for a 
total spacing of 3.0 (i.e. the spacing between the first 
and last source pixels) . 

20 

Turning now to Figure 2, a functional block diagram 
is shown according to a first embodiment of the invention 
for generating a stream of warped output pixels from a 
stream of source pixels in accordance with dynamically 

2 5 varying target increments. 

Firstly, the target increment (Tarlnc) is mapped via 
filter look-up table 21 to a specific FIR filter size as 
represented by a filter factor value based on the mapping 
30 shown in Table A, below. The FIR coefficients can be 

stored in a read-only table (i.e. ROM) or can be supplied 
by the user. Programming the coefficients gives the user 
the flexibility to tune the FIR filters. The filter size 
is chosen by indexing the filter look-up table based on 

3 5 Tarlnc. It is important to note that Table A is only an 
example of one implementation. This table could be 
extended to support any size of Tarlnc and generate any 
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filter size and filter factor. 



Table A - FIR Filter Table and Filter Factor 



Tarlnc (t) 


0<t<l 


Kt<2 


2<t<4 


4<t<8 


8<t<16 


16<t<3 


Filter Size 


3 Tap 


5 Tap 


9 Tap 


17 Tap 


33 Tap 


65 Ta 


Filter Factor 


0 


1 


2 


3 


4 


5 



Once Tarlnc is mapped, the largest filter on either 
10 side of a target pixel is chosen as the filter for that 

pixel. Figure 3 illustrates this step. The fifth line 
of Figure 3 (Tarlnc • Filter Mapping) shows the Tarlnc 
values from the fourth line (Target Increment) mapped to 
various FIR filters based on the mapping of Table A, 
15 above. This is represented schematically in Figure 3 by 

placing the largest filter size (i.e. number of taps) of 
the filter sizes appearing on each side of the target 
pixel inside the target pixel, indicating that it is the 
chosen filter for that pixel. In Figure 3, the third 
20 line (Target) shows the first target pixel being 

generated by a 9 Tap filter, the second target pixel 
being generated by a 9 Tap filter, the third target pixel 
being generated by a 3 Tap filter, etc. 



2 5 As discussed above, Tarlnc is mapped to an FIR 

filter size according to the mapping Table A, via filter 
look-up table 21, the output of which (i.e. filter factor 
value) is stored successively in A and B registers 22 and 
23, respectively, on successive clock cycles, thereby 

3 0 forming a pipeline. Comparator 24 determines which of 

two adjacent filters on each side of a target pixel is 
the larger, and stores a representation on this filter 
(i.e. the associated filter factor) in the left register 
25 connected to the output of comparator 24. The left 
35 kernel MAC 26 (i.e. multiplier accumulator) reads the 
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filter factor value stored in left register 25 and 
generates half of the target pixel value using the left 
half of the designated filter kernel. On the next clock 
cycle, the right kernel MAC 27 generates the rest of the 
5 target pixel value using the right half of the filter 

kernel and the centre tap stored in right register 28, 
and sums the resulting value with the output value 
received from left kernel MAC 26, resulting in the warped 
output pixel. 

10 

Before dynamic filtering can take place as discussed 
above in connection with comparator 24, the left and 
right registers 25 and 28 , and left and right kernel MACs 
26 and 27, the correct number of upsampled pixels must be 
15 generated. Referring again to Figure 2, the number of 

upsampled pixels is determined by the largest filter size 
stored before and after A and B registers 22 and 23. 

Comparator 29 determines the largest filter size 
20 (i.e. largest filter factor output from filter look-up 

table 21) and transmits that value to an upsample 
increment generator 30. The upsample increment generator 
30 receives and right-shifts the Tarlnc value by a 
predetermined number of bits equivalent to the received 
25 filter factor to produce an upsample increment value. 

For the example of Figure 3, the first target increment 
is 2.25, resulting in a filter factor of two which 
corresponds to a "divide-by-four" operation (i.e. right- 
shift Tarlnc by two bits resulting in an upsample 
30 increment of 0.5625). The upsample increment value (i.e. 

right-shifted Tarlnc value) represents the correct equal 
spacing of upsampled pixels for application to the 
required size of FIR filter (e.g. 9-Tap filter for Tarlnc 
= 2.25) . 

35 



Interpolator 31 reads the upsample increment value 
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and interpolates upsampled pixels from the source pixels 
on a continuous basis and with the spacing defined by the 
upsarople increment value, until an internal incrementor 
or counter value reaches Tarlnc. For the example 
discussed above and shown in Figure 3 , four upsampled 
pixels are generated with the spacing of 0.5625 relative 
to the source pixel spacing, for application to the 9 Tap 
FIR filter- The interpolator 31 may be implemented as a 
well known linear interpolator, a quadratic interpolator 
as disclosed in applicant's co-pending application serial 
number 08/172,065, or any other suitable design of 
interpolator. 

The upsampled pixels generated by interpolator 31 
are transmitted to the left and right kernel MACs 26 and 
27 which contain control logic to ensure that the correct 
coefficient is multiplied by the correct upsampled pixel. 
Thus, in Figure 3, line 2 (Upsample) depicts the largest 
possible number of required upsampled pixels for the 
different sized FIR filters implemented to generate the 
target pixels. The FIR filters shown operating between 
lines 2 and 3 of Figure 3 are implemented by left and 
right half kernel MACs 26 and 27 which multiply the 
correct coefficients by the correct upsampled pixel 
values. The internal control logic of the left and right 
kernel MACs 2 6 and 27 ensures that predetermined ones of 
the upsampled pixels are skipped for correct matching of 
coefficients and upsampled pixels. For example, whereas 
the 9 Tap FIR filter utilized to generate the first 
intermediate target pixel utilizes all of the upsampled 
pixels from line 2, the adjacent 3 Tap filter utilizes 
only the fifth, ninth and eleventh generated upsampled 
pixels, and skips the sixth, seventh, eighth and tenth 
upsampled pixels, etc. 

The control logic of left and right kernel MACs 26 
and 27 also generates the first target pixel on start up 
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of the system by doubling the right half kernel values 
for the first target pixel, as described in our copending 
Application No. 08/126,388. The last target pixel is 
generated by doubling the left kernel values. This 
5 doubling is required to ensure that all FIR coefficients 

maintain unity gain. 



Turning now to Figure 4, a functional block diagram 
is shown according to the preferred embodiment for 
10 generating a stream of warped output pixels from a stream 
of source pixels in accordance with dynamically varying 
target increments. 



Firstly, the target increment (Tarlnc) is mapped via 
15 filter look-up table 41 to a specific FIR filter size as 
represented by a filter factor value based on the mapping 
shown in Table A, above. The FIR coefficients can be 
stored in a read-only table (i.e. ROM) or supplied by the 
user. The filter size is chosen by indexing the 
20 coefficient table based on Tarlnc. It is important to 
note that Table A is only an example of an 
implementation. This table could be extended to support 
any size Tarlnc and generate any Filter size and Filter 
factor. 

25 

Once Tarlnc is mapped, the filter size is held in a 
register 42 and fed to the left and right kernel 
multiplier-accumulators 44 and 45. Figure 5 illustrates 
this step. The fifth line of Figure 5 (Tarlnc ® Filter 

30 Mapping) shows the Tarlnc values from the fourth line 

(Target Increment) mapped to various FIR filters based on 
the mapping of Table A, above. This is represented 
schematically in Figure 5 by placing the size of the 
respective half kernel of the assymetrical filter on each 

35 side of the target pixel, indicating respective half 

kernels for that pixel. In Figure 5, the third line 
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(Target) shows the first target pixel being generated by 
a 9 Tap filter, the second target pixel being generated 
by an assymetrical filter being one-half of a 9 Tap 
filter and one-half of a 3 Tap filter, the third target 
5 pixel being generated by a 3 Tap filter, and the fourth 

target pixel being generated by an assymetrical filter 
being one-half of a 3 Tap filter and one-half of a 5 Tap 
filter. 

10 As discussed above, Tarlnc is mapped to an FIR 

filter size according to the mapping Table A, via filter 
look-up table 41, the output of which (i.e. filter factor 
value) is stored successively in register 42. The left 
kernel MAC 44 (i.e. multiplier accumulator) reads the 

15 filter factor value stored in register 42 and generates 

half of the target pixel value and a first contribution 
to the center tap product value using the left half of 
the designated filter kernel and the center tap for the 
left-side half kernel. On the same clock cycle, the 

20 right kernel MAC 45 generates the other half of the 

target pixel value and a second contribution to the 
center tap product value using the right half of the 
filter kernel and the center tap for right-side half 
kernel. The right kernel MAC 45 then sums the first and 

25 second contributions to the center tap _ multiplication 
product and divides that summed value by two for 
obtaining an averaged center tap multiplication product, 
and finally sums the averaged center tap multiplication 
product with the left half of the target pixel value 

30 received from left kernel MAC 44 and the right half of 
the target pixel value, resulting in the final warped 
output pixel. 

As an alternative to summing the first and second 
3 5 contributions to the center tap multiplication product 

and dividing that summed value by two, the center tap 
coefficient for each half kernel may simply be halved 
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prior to storage in the associated coefficient ROM of 
left and right kernel MACs 44 and 45, or via the 
appropriate code where the invention is implemented in 
software, 

5 

The upsample increment generator 4 6 receives and 
right-shifts the Tarlnc value by a predetermined number 
of bits equivalent to the received filter factor to 
produce an upsample increment value. For the example of 

10 Figure 5, the first target increment is 2.25, resulting 

in a filter factor of two which corresponds to a "divide- 
by-four M operation (i.e. right-shift Tarlnc by two bits 
resulting in an upsample increment of 0.5625). The 
upsample increment value (i.e. right-shifted Tarlnc 

15 value) represents the correct equal spacing of upsampled 

pixels for application to the required size of FIR filter 
(e.g. 9-Tap filter for Tarlnc = 2.25). 

Interpolator 47 reads the upsample increment value 
20 and interpolates upsampled pixels from the source pixels 
on a continuous basis and with the spacing defined by the 
upsample increment value, until an internal incrementor 
or counter value reaches Tarlnc. For the example 
discussed above and shown in Figure 5, four upsampled 
25 pixels are generated with the spacing of 0.562 5 relative 
to the source pixel spacing, for application to the 9 Tap 
FIR filter. The interpolator 3 may be implemented as a 
well known linear interpolator, a quadratic interpolator 
as disclosed in applicant's co-pending application serial 
30 number 08/172,065, or any other suitable design of 

interpolator. 

The upsampled pixels generated by interpolator 47 
are transmitted to the left and right kernel MACs 44 and 
35 4 5 which contain control logic to ensure that the correct 

coefficient is multiplied by the correct upsampled pixel. 
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The FIR filters shown operating between lines 2 and 3 of 
Figure 5 are implemented by left and right half kernel 
MACs 44 and 45 which multiply the correct coefficients by 
the correct upsampled pixel values. 

5 

The control logic of left and right kernel MACs 44 
and 45 also generates the first target pixel on start up 
of the system by doubling the right half kernel values 
for the first target pixel, in the manner disclosed in 
10 Application No, 08/126,388- The last target pixel is 
generated by doubling the left kernel values. This 
doubling is required to ensure that all FIR filters 
maintain unity gain. 



15 Other embodiments and variations of the invention 

are possible. For example, according to a further 
alternative embodiment, the invention is implemented in 
software using C programming language. This alternative 
software embodiment is fully disclosed in Appendix "A" to 

20 this disclosure. Although the warping performed in 

accordance with the preferred embodiment has been 
described as providing proper Nyquist bandlimiting, 
special effect warping can also be implemented using the 
principles of the present invention without Nyquist 

25 bandlimiting (eg. soft focus warping). In other words, 
the pre-f iltering and interpolations steps can be 
implemented without Nyquist band-limiting. All such 
alternative embodiments and variations are believed to be 
within the sphere and scope of the claims appended 

30 hereto. 
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APPENDIX A 

17/ 

27/ Vertical and Horizontal Genesis Warping Software 1.0 

37/ 

47/This program warpc an input Targa file in both dimensions based on two 

57/ iooooyanymg vertical and horizonUl Warp Conuinmcot Files. 

67/ 

77/ Program Name: WARP.CPP 
87/ 

97/ Author James D. E. Gocl 
107/ 

1 17/ - Initial release of the software. 

127/ If you have any problems, please contact me. 

137/ 

147/ Copyright (c) ! 993 Genesis Microchip Inc. 

1 yj/ 200 Town Centre Blvd., Suite 400 

\6>jl Markham, Ontario 

17-7/ L3R8G5 

187/ (905)470-2742 

197/ 

20:// All rights reserved 

21:// 

227/ Modifications: 

237/ 

247/ Aug. 18, 1994 - James Goel 

23-y/ - Initial release 

267/ 

277/ Setup standard C definitions 
287/ 

29:0 include <stdio.h> 
30:#includc <stdlib h> 
3l:#inciude <math.h> 
32:*include <string,h> 
33:#indude <asserth> 
34:#indude <afich> 
35:#include "imagc.1T 
36:#indude "logiair 
37: 
387/ 

397/ Create a number of useful type definitions 
407/ 

4 Irtypcdef unsigned char BYTE; 

42 rtypedcf unsigned char PIXEL; 

43 rtypedef unsigned long ULONG, 

44 typedef unsigned tnt UINT; 
45: 

467/ 

477/ This header structure is used to define a 
487/ simple binary format for reading and storing 
497/ image data. 
507/ 

51 rtypedef struct image hdr 
32:{ 

53: // 

54: // Image Width and Height 

55: // 

56: unsigned intWidOu 
57: unsigned tnt Height; 
58: // 

59: // Image type 

60: // There are various Qiffcrenl Genesis define 

61: // images types: 

62: // 0 - Rod Pixels 

63: // 1 - Green Pixels 

64: // 2 -Blue Pixels 

65: // 3 -Luma Pixels 

66: // 4 - Chroma (U/V time multiplexed) 

67: // 5 - Black & White Pixels 

68: // 
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69: inl imagctype; 
70:} Genesis Header, 
71: 
727/ 

TMI Thia header structure is uaed to provide a 
747/ simple and oonaiatcnl method for parsing 
75:// command line parameters. 

767/ 

77:typcdef struct cmd line tag 
78:< 



79: char pchVTarlnc(13]7/ VcrticaJ warp containment field (VTarlnc) 

80: char pchHTarlncf 13 \JI Horizontal warp containment field (HTarlnc) 

81: char pchSrc[l3]; // Image to be warped 

82: char pchWarp[I3J; // Output warped image 



83:} tCmdLine; 
84: 
857/ 

86:// •••••*•••— ••••• 

877/ usage - Prints appropriate command line usage 

887/ 

897/ 

90: void UsageO 
91:{ 



92: // 

93: // Print out the proper usage 

94: // 

95: printf ( "ERROR: Incorrect command line parameters. VT); 

96: printf ( Usage: WARP -v [VcrticaJ Warp Oxuainrncnl File) "); 

97: printf ( "-h (Horizontal Warp ContainrnerU File]\n"); 

98: printf (" -w [Image to be Warped] 

99: printf ( "-o [Output Warped Filename JV); 

100: printf ( "Please do not use any wild card charactersAnW'); 



101:} 
1027/ 

1037/ •••••• 

1047/ ParseCmd - Parses the command line parameters 
1057/ •••••••••••••••••••••••••••••••••••••••••••••••• 

1067/ 

1 07: void Parse_Cmd( int argc, char •argvf], tCmdLine *ptCmd_Une) 
108: ( 



1 09: CFileExocption* the Exception = new CFiicExeeption; 

1 10: mt CurrArg-0; // Command line argument counter 

111: // 

1 12: // IF there arc no arguments supplied to the main line program. THEN 
113: // 

114: if (argc— 1) 
115: { 

116: // 

1 17: // Print usage information 

118: // 

119: UtageO; 

120: THROW( IhcException ); 

121: // 

122: //END IF -file arguments 
123: // 
124: } 
125: // 

126: // WHILE reeding through the command line arguments, 
127: // 

128: while ( ++CurrArg < argc ) 
129: { 

130: // 

131: //IF the input switches do not use dashes (-), THEN 

132: // 

133: if ( argv(CurrArgl(0] != V ) 

134: { 

135: // 

136: // Pring usage information 

137: // 

138: UsageO; 

139: THROW(theException); 
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140: // 

141: // END IF - input iwilcha 

142: // 

143: } 

144: // 

145: //SVVTTCHON:oofnrnAiklluv:«witch« 

146: // 

147: switch ( argv[CurrArg][l] ) 

148: { 

149: // 

130 // ON CASE: Vertical Warp Containment filename 

15li // 
132: caaeV: 
133: caaeV: 
154: // 

133- // Copy the warn containment filename 

136: " 

137 *rcpy< p<Cmd_Linc->pchVTarlnc, argv[++CurrArg] ); 

ISt: '•' 

139 it END CASE - Copy Straight Grid Filename 

160 . ' 

161: 

162: 

1$3 : /. OK CASE: Horizontal Warp Containment filename 

164: " 

165: «~ir: 

166: «a«V: 

167: // 

15g // Copv the warp containment filename 

> 6 * " ^ 

170: strcpy( psCmd_Uiie->pchHTarlnc t argv(++CurrArgl ); 

171: ' 

172: END CASE - Copy Straight Grid Filename 

173: " 

174: *«ak; 

175: " 

176: // ON CASE: Image to be warped 

177: /> 

178: oaae*W: 

179: caaeV: 

180: // 

HI- // Copy source image filename 

18* // 

1 g3 : strcpy< r*Crod_Unc->pchSrc, axgvf ++Curr Arg) ); 

184: A' 

185: // END CASE - Copy source image filename 

186: 7/ 

187: break; 

188: // 

189: // ON CASE: Output warped image 

190: // 

191: caacC: 

192* caac V: 

193: // 

194: // Copy output warped image filename 

195? // 

gtrcpy( psCmd.Line->pcfaWarp, argv[++CurrArgI ); 

197: // 

198: // END CASE - Output warped image 

199: // 
200: break; 
201: // 

202: // ON CASE DEFAULT: Show proper command line usage 

203: // 

204: default: 

205: // 

206: // Print proper usage 

207: // 

208: UsageO". 

209: THROW(theExoeption); 

210: // 
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//END CASE -default 

// 

break; 



If 

If END SWITCH 

II 
} 



II 

//END WHILE 

II 
> 



211: 
212: 
213: 
214: 
215: 
216: 
217: 
218: 
219; 
220: 
221: 
222:// 

213 J I End Pane Command Line 

2247/ 

225:) 

2267/ 

2277/ • 
2287/ Main program loop 

2297/ ••••••••••••••••••••••••••• • • 

2307/ 

23 1 :void main(int argc, char •argvfj) 
232:{ 

// Exception object for debugging 
CFileException* thcExccption - new CFileException; 



233 
234: 
235: 
236: 
237: 
238: 
239: 
240: 
241: 
242: 
243: 
244: 
245: 
246: 
247: 
248: 
249: 
250: 
251: 
252: 
253: 
254: 
image 
255: 
256: 
257: 
258: 
259: 
260: 
261: 
262: 
263: 
264: 
265: 
266: 
267: 
268: 
269: 
270: 
271: 
272: 
273: 
274: 
275: 
276: 
277: 
278: 
279: 
280: 



Image Srclmage; 
Image Tar Image; 

Pixel VUpPixel; 
Pixel VTihPix; 

LeftKemel VLeAKemc!; 

RightKemel VRightKemel; 
FILE •pfVTarlnc; 
UINT iVTarlncWidth; 
UINT iVTarlncHeight; 
UINT iTarTypc; 
float •pVTarlncLinc; 
float •pVTarlncLineTmp; 

tCmd Line aCmd Line ; 

LOOKUP VFihLookup; 

IncGen VUpInoGen; 

Interpolate Vlnterpolator, 
BYTE •pSroCol; 
BYTE •pTarCol; 
BYTE •pTarColTmp; 
UINT iColumn; 

i Count, 



// Source image to be warped 

// Output warped image 

// Vertical Up sampled pixel 

// Vertical Filtered pixel 
// Vertical Left FIR Filter Kernel MAC 
// Vertical Right FIR Fitter Kernel MAC 

// File pointer to the Vertical Tarlnc 

// Width of the the containment file 
// Height of the containment file 

// •*• Reserved for future expansion ••• 

// Pointer to a vertical warp Tarlnc 
// Copy of the previous pointer 

// Holds the command line parameters 
// Maps the Vert Tarlnc to FIR fitter size 

// Object to generate a Vertical Up Sample increment size 
// Object to generate an interpolated vertical up-sampled pixel 

// Pointer to a source image column 

// Pointer to a target image column 

// Copy of the previous pointer 

// Counts the number of columns in the input 



UINT 

// 

// Print introduction 

// 

Iprintf ( stderr, "\n\nWARP FAST! Generates a warped Targa file. Aug-23-1994W); 
fprintf ( atderr. *Copyrighl<c) 1994 Genesis Microchip. All Rights Reserved^"); 
fynatf ( stderr. "Author. James GoeI\n\n"); 

// 

// Parse the command line parameters 
// Load the command line parameters into the command line 
//a 

// 

TRY 



// Counts the number of output pixels generated 



Pane_Cmd( argc, argv, &*CmdJLine ), 



( 
) 

CATCH( C Exception, theException ) 
{ 

TRACE( "Exception: Error during command parsing phase. Jn" 

«*(!); 

) 

END CATCH 

// 

// •••••.••«••••••,•»«••,.. 

// VERTICAL WARPING ALGORITHM 
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282: // Open the Vertical Warp Containment file (Tarlnc) file 

283: // IF the Tarlnc file could not be opened, THEN: 

111: K (pfVTarlnc-fopen( ^od_Line pcfa VTarlnc, "rb" )) « NULL ) 

286: { 

287: // 

288; // Print an error menage and abort 

2^' f fpnrtf{ atderr. 'Could not open the Vertical Warp Containment CleAn"); 

291: fyfiitti{ stderr. "Please correct this problem and try again.\n"); 

292: ewt(0); 

293: // 

294: //END IF 

295: // 

296: ) 

297: // 

298: // Read the Width and Height of the Vertical Tarlnc file 

299! // 

300: 6«ad( AiVTarlnc Width. sizeof^UINT). 1, pfVTarlnc ); 

30 1 : m=ad( &i VTarlncHcight, sixeof^UINT). 1, pfVTarlnc ); 

302: // 

303: //Read the type of Vertical Tarlnc file 

304: // (Used for future expansion) 

305: // 

306: fread( AiTarTypc. iizeof(UINT). L pfVTarlnc ); 

307: // 

308: // Allocate memory to hold one line of Vertical Tarlnc values 

309: //IF the line could not be allocated. THEN: 

3 ll! if( (pVTarlncUne - (float «>nAlk>o<«i2iM^no*t)^VT*rIncWidlh)) — NULL) 

312: { 

313: // 

314: // Print error message and abort 

315: // 

316: fynx*f{ Aden, "Could not allocate memory for the Vertical Tarlnc lincW); 

317: odt(l); 

318: }; 

319: // 

320: //Open and load the input Targa File 

321: // 

322: TRY 

323: { 

324: Srclmage.Load( sCmd_Line pohSrc ); 

325: ) 

326: CATCH( CExocptioo, thcExceptiou ) 

327: { 

328: TRACE( "Error while loading the Source Image.ta" ); 

329: fdomealK): 

330: exh(l); 

331: ) 

332: END CATCH 

333: // 

334: // Verify that the width of the Coouinmcnt file n 1 less 

335: // than the width of the aource Targa file. Thii enure 

336: // that only the spaces betwoen the iriput Targa lines get 

337: //resized. 

338: // 

339: aasert( iVTarlncWidlh — Srcimage.HeightO - 1 X 

340: // 

341: // Verify that both VTarlnc and the Source image have the 

342: // same number of columns 

343: // 

344: asscrt( i VTarlncHeight — SrcImage.WidthO X 

345: // 

346: // Print sUtuj meaasgc 

347: // 

348: fprirril atderr. "Loading source image.. ltt*]\n\ sCmd Linc.pohSrc ); 

349: // 

350: // Allocate enough memory for the output Vertical Wsrp file 

351: // 
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J J A. 


TRY 




J JJ. 


{ 








Tarlmage.Crcate< Srclmage. WiduSQ. Srclmage.HrightO ). 


J J J. 


) 




J JO. 


CATCH( CExDCfXioa, theException ) 


J 57. 


( 




J3B. 




TRACE( "Error while allocating space for the Target Warp ImageAn" ); 


359: 




fcW*U0 ( 






«□«.); 


JO 1 . 


) 




JO*. 


END CATCH 


J6J. 


If 




364: 


// LOOP for the number of vertical column* in (he Containment Field 


Jo J. 


It 




vex 


fbr( iCc4umn=0; tColumn < iVTartncHeight; iCoIumn++) 


JO ' . 


( 




JoJ 




If 


1XQ 
JOT 




// Read in a Vertical Tar Inc Line 


J *v 




/' 


J ' 1 




fr«a4 pVTar Incline, tizoo^ftoat), iVTarlnc Width, pfVTarlnc ); 


1T2 






J * 1 




'■ Ummd a souroe image column 


I'M 










pSmiCat - SrcImage.GetColumn(iCoiumn); 


J 7% 






JT7 




I *mA fee source image column into the Vlnicrpolator 


j /• 










« wMrpoutor.>eC>rc( porcCol ), 


310 






311 




' irti a target tmage column that will be loaded with 


js.* 




taanwi tfnAf*^ filtrate 

'■w^*o image pixcu. 


383 










• TwOal • TArlmaiye GelCnliimnfiCnlumnY 

^ * ■ «i III Wgw.VMWlllllll^lVvllUIUI J 9 


385: 






386: 




f bwuauxe the Vlnterpolator and make copies of the 


387: 




// ^oatm to both the Tar Inc and Target image files 


388: 




a 


389: 




Vlflurpolator.ReaetO; 


ion 




Vliacrpolalor.iTarlncWidth - iVTarlnc Width; 


101 - 




Y * 1 •» UMi4iWi A U(|S VTA 4U MWa sksMPSj 


392: 




pTarCotTmp - pTarCol; 


393: 




// 


394: 




H rYrt working line number 


395: 




// 


396: 




tprtndt stderr, "Working on column.. .{%dpo". iColumn ); 


397: 




// 


398: 




//Reaet the Right Vertical MAC 


399: 




// 


400: 




VTRightKemcLReaetO; 


401: 




// 


402: 




// LOOP for the number of Tar Inc values in the Vertical Conuinment Field 


403: 




// 


404: 




for( iCount-O; iCount < iVTarlnc Width; rCount++ ) 


405: 




( 


406: 




// 


407: 




// Reaet the Left Vertical MAC 


408: 




// 


409: 




VUftJCerneLReaetO; 


410: 




// 


411: 




// Calculate the correct Vertical Fitterfactor based on Tar Inc 


412: 




// 


413: 




VFihLoofcup fTarlnc - •pVTarlncLineTmp; 


414: 




// 


415: 




// Load the current VTarlnc into the Vertical Upaample Increment Generator 


416: 




// Load the correct Vertical Fitter Factor into the Uptamptc Increment Generator 


417: 




// 


418: 




VUpInoGaufTarlnc - •pVTarlncLineTmp; 


419: 




VX^plncOcfLFiltFact - VFihLcokup.OutO; 


420: 




// 


421: 




// Load the current Tar Inc into the VlnterpoUior 


422: 




// Load (he Filter Factor into the Vlnterpolator 
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423. // Load the Upsample Increment into the Vlnterpotator 

424: " 

425 : VInterpoiaioriTarinc - ^VTarlncLtneTmp++; 

426: Vlnterpolator.FUlFact - VFillLookup.OutO; 

427 VlaLapoUior.UpInc - VUpIncGai-OutO; 

428^ // 

429 // Load the Finer Factor into the MAC* 

// 

431 VLeftKeroeLFUtFaot - VFihLookup.OutO; 

432 VRightKeroeLFiKFact - VFihLooicup.OutO. 
433: // 

434- // 0 7 the output pixel is a top edge pixel, THEN: 

43* // 

436: "T iCount— O > 

437: // 

43g. // Double the Left filter kernel coefficients 

43* " 

440: VRightKemclJ>ouWe»TRUE; 

441: // 

442 : // ELSE IF the output pixel is a bottom edge pixel, THEN: 

443! // 

444 ebe if{ iCounl— iVTarlncWidth-l ) 
445: // 

446 // Double the Right fiher kerne! coefficients 

447: // 

448: VLeftKemel.Double«TRUE; 

449: " 

450: // ELSE if the output pixel is somewhere in the middle of the image. THEN: 

451: // 

452: else 

453: { 

454: // 

455. // Do not double either Right or Left kernel 

45* // 

457: VLeftKcmel.Double-FALSE; 
458: VRightKemeKDouble-FALSE; 

459: } 

460: // 

46L // END IF - output pixel is on an edge 

462: // 

463: // 

464 : // Calculate the center tap Kernel values 

465^ // 

466: VLeflKemeI.KMC(VInlerpoUtor.UpPix); 

467: VRightKcmeI.KlAC( VlnterpoUtor.UpPix ); 

468: // 

469: // LOOP through the other FIR filter taps 

470: // 

471: while< VlnterpoUtor.CenterTapO f-TRUE ) 

472: { 

473: // 

474: // Generate an interpolated pixel 

475- // Filter the interpolated pixel 

476: " 

477: VUpPixel - V]i*efpc4ator.Cknerate<); 

478- VLeflKemeLMAC( VUpPixel ); 

479^ VRighlKemel.MAC( VUpPixel ); 

480: // 

481: // END LOOP - other FIR fitter taps 

482: // 

483: ): 

4*4: // 

4g5: // Write the Right Kernel Output Pixel to the output line store 

486: // 

487: VFihPix - VRighlKemeLFiherPixelO; 

488: •( P TarCoITmp+^) - VFihPix. Red; 

489: •(pTarCcJTmjrM.) - VFUtPixGm; 

490: •(pTarCodmjrM') - VFihPix. Blu; 

491: // 

492: // Check to make mrc that the Target Unc stays within output image range 

493! // 
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* , ; : 



45M: 






aO V 




// 


AQJZ- 






A<il- 

i . 




// 






VRjfftiiKtfrud AflcmnRfl<fr»VL^ftK«My AocmnRed: 






VP iofrf K<wr^l AmrroOrr*" VT ^fHCf*r^l AraimTrm' 






VRJehtKariei AccumBlu«VL£ftKa«l AccumBlu* 


501: 






302. 




// cinl* L*u«jr — i ix loc vmiucs crorn waauuaiicfa rue 


303: 










* 


305: 






306: 




// write Lcn Kernel vjuipui nxel to toe output une norc 


307: 






(AO, 

308: 




vrinrix ■ vijeiuAxmeLruxerrtxei(j, 


309: 




(P I trCol imp* * J ■ vriiirULKea, 


sin- 




</ n T. f ^V t |T mn i »\ « VFihPir On' 
1 waii i mp ' * / * r nix IA.VJTU, 


311. 




*/nT»rTiilTmn i * \ m WthPiv Rln- 

y> i u\«oi i mpTTj ■* vnnruLinu, 


JU. 




// 


(11, 

31 J. 




// Check to make sure that the Target Line stays wvthm range 






// 


J I J . 




men^ ^mii i / \jf i aiwot i ni|>p i «r\«oi* i ) i ■rurugc. w lauw^j j 


JlO, 




// 


517: 




// wnie me vjuuxn rtxci une to otaic 


518: 




// 


519: 




*T~D V" 


520: 




{ 


521: 




TarImage.PutCoiuinn( pTarCoU iColumn ); 


522: 




} 


523: 




CATCH( CException. theException ) 


524: 




( 


525: 




TRACER "Error while writing • warped output image line.\n M ); 


526: 




fcloseallO; 


327. 




exrt(l); 


32B. 




) 


327. 




END CATCH 


530: 


If 




531: 


If END LOOP - vertical columns 


532: 


II 




533: 


) 




534: 


II 




535: 


II Save vertically warped Targa file 


536: 


II 




317- 
j j / . 


TRY 




J JO. 


{ 




jjy. 




TarImage.Save< "verttga" ); 


540: 


) 






CATCH( CException, theException ) 




( 








TRACE( "Error while saving the vertically warped tmageVn" ); 


544: 




fcloaeallO; 


545: 




exrt(l); 


546: 


> 




547: 


END CATCH 


548: 


II 




549: 


II Cloae the Vertical Tar Inc file 


550: 


II 




551: 


fcloae<pfvTaflnc). 


552: 


II 




553: 


//Print* 


atua message 


554: 


II 




555: 


fyrintfl stderr, "Finished Vertical Warping.. An"); 


556: 


fprintf^ stderr, "Starting Horizontal Warping. Jo"); 


557: 


II 




558: 






559: 


II Horizontal WARPING ALGORITHM 


560: 


If •••••• 


•*•••••••«•••••••••••••••••«••••••«*•**••••*«••«•••••*••«*•••••• 


561: 


II 




562: 


II Define all Horizontal warping variables 


563: 


II 




564: 


Pixel 


KUpPtxel; //Horizontal 
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565: Pixel 
566: Image 
Horizontally warped 
567: Image 

LefUCeroel 
RightKeroet 
FILE 
UtSJ 
UINT 
float 

float 

LOOKUP 
InoGco 
Interpolate 
BYTE 
BYTE 
BYTE 
UINT 
l image 
// 



HfihPix; 
HSrclmage; 

HTarlmage; 
HLeflKernel; 
HRighlKcroel; 

•pfHTarlnc; 

tHTarlncWidux 

iHTarlncHeight; 

•pHTarlncLine; 



// Horizontal Filtered pixel 

// Vertically Warped Source image to be 

// Output warped image 
// Horizontal Left FIR Fitter Kernel MAC 
// Horizontal Right FIR Filler Kernel MAC 

// Pointer to Horizontal Warp Tarlnc 
// Width of the the corUaioment file 
// Height of the containment file 

// Pointer to a specific Horizontal warp containment (Tarlnc) 



•pHTarlncUneTmp; //Copy of the previous pointer 



HFiltLookup; 
HUpIncGcn; 
HlmcrpoUtor; 

•pSrcRow; 

•pTarRow. 

•pTarRowTmp; 

iRow; 



// Object to map the HorizonUl Tarlnc to the correct FIR filler size 

// Object to generate a Horizontal Up Sample increment lize 
// Object to generate an interpolated HorizonUl up-aampled pixel 
// Pointer to a source image row 
// Pointer to a target image row 
// Copy of the previous pointer 

// Counts the number of rows in the 



// Open the Horizontal Warp Containment file (Tarlnc) file 
//IF the Tarlnc file could not be opened, THEN: 

// 

ifl (pfHTaruic=fopcn( sCmd Line.pchHTarlnc. "m" )) = NULL ) 
( 

// 

// Print an error message and abort 

// 

fprintfl stderr, "Could not open the HorizonUl Warp Containment file.\n"); 

fixing stderr, "Please correct this problem and try again.\n"); 

cxh<0); 

// 

//END IF 

// 

> 

// 

// Read the Width and Height of the Horizontal Tarlnc file 

// 

fread( AiHTarlncWidth, sizeof(UINT), 1. pfHTarlnc ); 
fiead( AiHTarlncHeight, rizeo^UINTX 1, pfHTarlnc ); 

// 

// Read the type of Tarlnc file 
// (Used for future expansion) 

// 

fread( AiTarTypc, sizeoffUINT), I, pfHTarlnc y t 

if 

// Allocate memory to bold one line of HorizonUl Tarlnc values 
//IF the line could not be allocated, THEN: 

// 

if( (pHTarlncLine - (float *>iiaJloc<aize<>((flcat)*iHTaruicWidm)) — NULL) 



// Print error menage and abort 

// 

fyr\ntf{ stderr, "Could not allocate memory for the HorizonUl Tarlnc linc.\n"); 
exitd); 



// 



// Open and load the Vertically warped Targa file 

// 

TRY 



HSrclmage. Load( "verUga" ); 



( 
) 

CATCH( CExcepuon, thcExccption ) 

TRACE( "Error while loading the Source image.Nn" ); 

fdoseallO; 

exhtl). 
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633: ENDCATCH 

634: It 

635: // Verify thai the width of the Conuinmrrt file u 1 leas 

636: // than the width of the source Targs file This ensures 

637: // that only the space* between the input Targa lino* get 

638: // raized 

639: // 

640: asaort< iHTarlncWtdlh — HSrcImagc.Width() - 1 ); 

641: // 

642: // Verify that both HTarlnc and the Souroc image have the 

643: // same number of Una 

644: // 

645: asaert( uTTsrlncHeight — HSrcLmage,HeightO X 

646: // 

647: // Allocate enough memory for the output Vertical Warp file 

648: // 

649: TRY 

650: { 

65 1: HTaruriage.Create< HSrcImage. WidthO, HSrclmagcHeightO ); 

652: } 

653: CATCH(CExcepUorwtheExcepuan) 

654: { 

655: TRACER "Error while allocating space for the Target Warp Image.V ); 

656: fclceealK); 
657: exit(l); 

658: } 

659: END_CATCH 

660: // 

661: // LOOP for the number of Horizontal rowi in the Containment Field 

662: // 

663: for( iRow-0; iRow < iHTarlncHcighl; iRow++) 

664: { 

665: // 

666: // Read in a Horizontal Tar Inc Line 

667: // 

668 : tread( pHTarlocLinc, iizeof(noat), tHTarlnc Width. pfHTarlnc ); 

669: // 

670: // Read a source image row 

671: // 

672: pSrcRow - HSrcImage. GctUxWi Row), 

673: // 

674: // Load the souroc image row into the H Interpolator 

675: // 

676: HinterpoUtor.SctSrc( pSrcRow ); 

677: // 

678: //Get a target image row that will be loaded with 

679: // warped image pixels. 

680: // 

68 1 : pTarRow - HTarlmage.GetUneftRow); 

682: // 

683: // Initialize the Hinterpolator and make copies of the 

684: // pointers to both the Tar Inc and Target image files 

685: // 
686: H Interpolator. RosetO; 

687: Hlnierpolator.iTarlnc Width - iHTarlnc Width; 

688: p HTarlnc LineTmp - pHTarlncLinc; 

689: pTarRowTmp - pTarRow, 

690: // 

691: // Print working line number 

692: // 

693: fprintfl stderr, "Working on row...[%d]\fl-, iRow ); 

694: // 

695: // Reset the Right Horizontal MAC 

696: // 

697: HRigrtiKernd.ResetO; 

698: // 

699: // LOOP for the number of Tarlnc values in the Horizontal Containment Field 

700: // 

701: for< iCount-0; iCount < iHTarlncWidth; iCount^+ ) 

702: { 

703: // 
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// React the Left Horizontal MAC 

// 

HLrftKcroeLRadO; 

// 

// Calculate the oocroct Horizontal Fihcrfactor baaed on Tarlnc 

// 

HFihLookupiTar Inc - •pHTwInclincTmp; 

// 

// Load the current HTarlac into the Horizontal Upaample Increment Generator 
//Lc^thecorrecxHorizcnlalFihe^^ 

// 

HUpIncGen/Tarlnc - •pHTtrlncUneTmp; 
HUplnoGcn.FihFACt - HFihLookup.OutO; 

// 

// Load the current Tarlnc into the HlrterpoUtor 

// Load the Filler Factor into the HIntcrpolalor 

// Load the Upaample Increment into the Hlnterpolator 

// 

HlnterpoIaloriTarlnc * ^HTirlncUneTmp^ ; 
Hlnterpolator.FihFact = HFiltLookup.OutO; 
HlnterpoUtor.Uplnc * HUplncGertOutO. 

// 

// Load the Filter Factor into the MAC* 

// 

HUftKerneLFihFact « HFihLookup.OutO; 
HRightKernel.FiKFact = HFiKLookup.OutO. 

// 

//IF the output pixel U a top edge pixel, THEN: 

// 

uTiCount— 0) 

// 

// Double the Left filter kernel coefficients . 

// 

HRighlKemcl.Doublc-TRUE; 

// 

// ELSE IF the output pixel it a bottom edge pixel, THEN: 

// 

clie if( iCount==iHTaf IncWidth-I ) 

// 

// Double the Right Cher kernel coefficient* 

// 

HLeftKeroel.Double=TRUE; 

// 

// ELSE if the output pixel it somewhere en the middle of the image, THEN: 

// 

eUe 

{ 

// 

// Do not double either Right or Left kernel 

// 

HUftKerneLDouble-FALSE; 
HRigtuKerncl.Doublc= FALSE, 

} 

// 

//END IF - output pixel u on an edge 

// 
// 

//Calculate the center tap Kernel values 

// 

HLcftKernel.MAC( HlnterpolatorUpPix ); 
HRightKemeLMAC( Hutterpoiator.UpPix ); 

// 

// LOOP through the other FIR filter tap* 

// 

*tiile< Hlrterpolatc<.C^ntcrTapO !- TRUE ) 
( 

// 

// Generate an interpolated pixel 
// Filter the interpolated pixel 

// 

HUpPixcl - HlnterpcUalor.GencratcO. 
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775 




HLeAKerneLMAC( HUp Pixel ); 


776 




HRighlKenieLK(AC( HUp Pixel y. 


777 




II 


778 




II END LOOP - other FIR fiher tap* 


779 




// 


780 




); 


781 




// 


782 




//Write the Right Kernel Output Pixel to the output line store 


783 




// 


784 




HFittPix - HRightKernetFiherPixelO; 


785 




•(pTarRowTmp-M- ) - HFiKPix.Red; 


786 




•(pTarRowTmp++) - HFihPix-Gm; 


787 




•(pTarRowTmp-H-) - HFihPix.Bhi; 


788 




// 


789 




// Check to nuke sure that the Target Line stays within output image range 


790 




// 


791 




aaaert( (UINT) (pTarRowTmp-pTarRow-1) < HTarImage.Width0*3 ); 


792 




// 


793 




//Transfer the Left Kernel to the Right Kernel 


794 




// 


793 




HRighlKewLAixumRed-HUftKerncLAocumRed; 


796 




HRightKefnel.AocumGm-Hl^ftKernel.AecumG^ 


797 




HRjghtKcrncl.AocumBlu-HUftKcrncl.AccumBlu, 


798 




// 


799 




// END LOOP - Tarlnc values from Containment File 


800 




// 


801 




} 


802 




// 


803 




// Write Left Kernel Output Pixel to the output line store 


804 




// 


805 




HFihPix - HLeftKemcl-FiherPixelO; 


806 




•(pTarRowTtnp++) - HFittPixRed; 


807 




•(pTarRowTmp-M-) - HFihPix.Gm; 


808 




*(pTarRowTmp++) - HFihPix. Biu; 


809 




// 


810 




// Check to make aurc that the Target Line slays within range 


811 




// 


S12 




aaaert( (UINT) (pTarRowTmp-pTarRow-|) < HTarlmage. WiAh0*3 ); 


813 




// 


814 




// Write the Output Pixel line to disk 


815 




// 


816 




TRY 


817 




( 


818 




HTar Image. PutUne< pTarRow, iRow \ 


819 




> 


820 




CATCH( CException, the Exception ) 


821 




{ 


t22 




TRACER "Error while writing a warped output image line.Va" ); 


823 




fcWaliO; 


824 




«at(l); 


825 




} 


826 




END CATCH 


827 


II 




828 


II END LOOP - horizontal lines 


829 


II 




830 


> 




831 


II 




832 


II Save warped Targa file 


833 


II 




834 


TRY 




835 


{ 




836 


HTarlmagc.Save< sCmd Une.pchWaip X 


837 


) 




838 


CATCH( CException, theExoeption ) 


839. 


{ 




840 




TRACE( "Error while saving the final Warp imageW \ 


&41. 




CcloaeaUO; 


842: 






843: 


} 




844 


END CATCH 


*45: 


II 
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BAA 


// Free all allocated memory 


a ji. 


11 


W: 


fr«< pVTar Incline \ 


849: 


frtx< pHTarlncUne ); 


850: 


SrelmagaFreeManoryO; 


851: 


Tarlmage.FreeMenioryO; 


852: 


HSccImage-FreeMemoryOi 


853: 


HTaflmagcFreeMemoryO; 


854: 


fcloaeallO*. 


855: 


// 


856: 


//PriiAitatottmeaaage 


857: 


// 


858: 


fpnr*f{ rtderr, "Program completed wece«irully!\n" y. 


859:// 




*6QM END - Warptn* FAST! 


UVJl 




862:} 




863: 
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Id l Image Clan - Enrapsnlilm the Leadfcoots Image Functions 



4//Thif image cUaa greatly improvce the LeadTooU user mtrrfanr by 
5JI providing a simplified and consistent ad of C++ member functions 
6'M for manipulating images. 

I'M 

%M Program Name: IMAGE.CPP 

9-y/ 

10:// Author James D. E. Gocl 
li'M 

12M - Initial releaae of the software. If you have any problems, please oontact roe. 
13:// 

XA'JI Copyright (c) 1993 Genesis Microchip Inc. 

\5M 200 Town Centre Blvd.. Suite 400 

16 j/ Markham, Ontario 

17M L3R8G5 

ISM (905)470-2742 

19M 

20 j/ All rights reserved. 
21M 

22:// Modifications: 

23:// 

7AM Aug 5. 1994 - James God 

25 M - Initial release 

26M 

27:#include-iroagc.h- 

28:*uiclude<stdlib.h> 

29:#includc <aasert_h> 

30:#indude<afiLh> 

31: 

32M 

337/ Error codes 
347/ 

35:WefineNAME_SIZE 1 
36:#defineSAVE_ERR 2 
37:#defineREAD ERR 3 
38:#define MEM ERR 4 



44M 

45 .-void Image:: Errcr( tnt crrnttm ) 
46:{ 

47: // 

4*: // SWITCH ON: Associated error 

49: // 

50: iwitcn<erroum) 




42M error - Prims appropriate error message and aborts gracefully. 




31: 
52: 
33: 
54: 
55: 
56: 
57; 
58: 
59: 
60: 
61: 
62: 
63: 
64: 
65: 
66: 
67: 
68: 
69: 
70: 
71: 



// 

case NAME SIZE: 

II 

II Prmt out a file name size error message 

// 



( 



// 

// END CASE: command line 

// 

break; 

// 

// ON CASE: Leamools could not cave (he output image 

// 

case SAVE ERR: 



// ON CASE: Command line parameter error 



// 



// 

// Print out save error message 

// 



fpriiHf ( stderr, "ERROR: One of the input file names exceeds the legal 8 J file name size.\n"); 
fprintf ( stderr, "Please reduce the size of the file name and try again.\n" ); 
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72: fronts aiderr. 'ERROR: Could ooc save the output image filc.\n"); 

73. £pnnlf( Hxkrr, Tbc disk nuy be full or corrupted "); 

74; fprintf( stderr. "Please correct these condition* and try agsin.\n-); 

73! // 

76: //END CASE: save error 

77: // 
78: break; 
79: // 

SO: // ON CASE: Could not read the input image 

ti: II 

12: ca»c READ ERR: 

83: // 

$4: // Print out load error menage 

ty // Stop execution 

96: II 

f7: frrinxf( stderr, "ERROR: CouW not read one of the the input filesAn"); 

gS; %nntf( stderr. "The disk nuy be corrupted"); 

89; frhnd{ stderr, "Please correct these condniom and try agaia\n"); 

90: " 

91: /' EVD C Alt load error 

92: '-■ 
93: arsaa, 

94: *■ 

95: /' C Kit Could not allocate the required memory 

96: 

97: «m> ut M IUL 

98: 

99: frwc out error 

100: 

101: f^«^ stderr. "ERROR: Could not allocate the required rnernory.Nn*); 

102: atdcrr. "Please free up aome DOS memory and try again.\n"); 

103: 

104: /'END CASE memory error 

105: * 

106: 

107: /' 

108: /' ON DF FAULT: A fatal internal error has been generated 

109: /' 
110: default 
111: // 

112: // Print out UiaJ mm*! error 
113: // Stop cxrouuon 
114: // 

113: pnrtffTERROR: A fcUal internal error has occurred. Please call James GoeUn"). 

116: // 

117: // END DEF AULT: fmtal error 

118: // 
119: break; 
110: // 

121: // END SWITCH: associated error 
122: // 
123: } 



124:// 

125:// END ERROR 

126// 

127:> 

128: 

129:// 

130-7/ Default Constructor 
1317/ 

132:Image::Image(votd) 
133:{ 

134: FlagAUocLine - FALSE; 
135: FlagAllocCol - FALSE; 
136// 

137// END Default Constructor 

138// 

139:} 

140: 

141// 

142// Default Destructor 
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161: delete pUnc; 



143:// 

1 44:Iinage::Mmage<vo<d) 
U5:{ 

146: TTRACE( "Freeing BitaM|)...\n"X 
147: L_FrecBitm«p( &Image Bitmap ); 
14*: delete pLine; 
149:// 

1507/ END Default Deatntctor 

151:// 

152:} 

153: 

154:// 

155:// Free Bitmap* 
156:// 

157:void Imigc::FrecMcmocy(void) 
158:( 

1 59: TRACE( "Freeing Bitmap and memory „>a m ); 
160: L_FreeBitmap( & Image Bitmap ); 



162: delete pCol; 
163:} 
164: 
165:// 

166:// Construct an empty image 
167:// •••••••••*•••••*••••••** 

168:void Image: :Create( U1NT ilmgWidth, UINT ilmgHeight ) 
169:( 

CFileExocption •the Exception ** new CFileExceptiocu 

// 

f Initialize the bitmap to new width and height 
//IF the bitmap wat HOT initialized auooeaaiuUy, THEN: 

// 

LJnhBrtmap( &lmage Bitmap, ilmgWidth. ilmgHeight, 24); 

// 

// Allocate memory 

//IF the memory cannot be allocated, THEN: 

// 

lf{ L_AllocaieBitroap( AlmageBitmap, TYPE_CONV ) 1- SUCCESS ) 
i 

// 

//Printi 

// 

Enor( MEMERR X 
L_FrocB*tmap( ftlmageBftmap ); 
delete pLine; 
THRO W( thcExoeption ); 

// 

// END IF - Allocate the bitmap memory 




) 

t Width - 
tHcigbt- 

// 

// Initialize aome baat c Bitmap propertiea 

// 

InugeBtonap.VtcwPerapeotive - BOTTOMLEFT; 
Image Bhmap.Order - ORDERBOR; 
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2,4: i( L_l^QAScreon< pcnlnput. *mageBitmap. OUTPUT_BITMAP, TYPE.CONV, 0. 0)!- SUCCESS ) 
215: ( 

216: // 

217- // Print error message and abort 

218^ // 
219: Error( READ_ERR); 

220: THROW(thcExoeptk»); 
221: // 

222: // END IF - load Targa image 
223: // 
224: } 
223: // 

226: //Copy the vndmaiidbriglBmf^^ 
227: // 

228: iWidth - ImageBttmap/Width; 
229: Weight -ImageBitmap.Hcight; 
230-7/ END - Load image file 
231:// 
232:} 

2337/ ••••••••*••—••••• 

2347/ Save an Image FUe 

2357/ 

236:void Image::Save< char *pchOitput ) 
237:{ 

238: CFUeExocption •theException - new CFilcEtocpUon; 
239: // 

240 : // Save Targa Image 

241: // IF the image could re* be uved, THEN: 

242: // 

243 : U( L_SavcTG ABitmap( pchOutput, AlmageBhrnap, 24)!- SUCCESS ) 
244: { 

245: // 

246: // Print error message and abort 

247: // 

248: Error< S AVE_ERR X 

249: L_FreeBitmap( &ImageBitmap ); 

250: delete pLme; 

251: THROW(theExoeptioa); 

252: // 

253: //END IF- Save Targa image 
254: // 
255: } 
2567/ 

2577/ END - Save image file 

2587/ 

259:) 

2607/ ••••••••••••••••••••••• 

2617/ Return the Image Width 
2627/ 

263:UINT Image::Width( void ) 
264:{ 

265: retum< iWidth); 
2667/ 

2677/ END - return image width 

2687/ 

269:} 

2707/ 

271 7/ Return the Image Height 
2727/ ••♦•••••••••••*•••••••• 

273:UINT Image:: Height void ) 
274:{ 

275: fotum< iHetght); 
2767/ 

2777/ END - return image Height 

2787/ 

279:} 

2807/ •••••••••••*•*••••*••••••• 

2817/ Get a line from the Image 
2827/ *•••••*•••**•••*••••*••**• 

283:BYTE •image::GcOinc< UINT iUneNum ) 
284:{ 
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285: CFdcExoeption •the Exception - new CFUcExeeptioo; 

286: «Meri< iUneNum < (UINT) ImageBitnup.Hcight ); 

287: // 

288: //IF * line has not been allocated, THEN: 

289: // 

290: if( FlagAUocUnc — FALSE) 

291: ( 

292: // 

293: //AUocalc the Line 

294: // 

293: pUne - new BYTE(IinAgeBitmap.Byte« Per Line], 

296: FlagAllocLine - TRUE; 

297: // 

298 // EKD IF - a line baa been allocated 

299. // 

300 ) 

301 ft 

302 it IF mm on^a m m the bottom left comer, THEN: 
301 /' 

304 * fanjiflaaia V^Pcnpoctive — BOTTOMLEFT ) 

30) ( 

)04 ' AAwat the row number based on the origin of the image 

307 U awHwa ■ Image Bitmap. Height - ilincNurn - 1; 

30t / 

309 / tSV If - an^a m bottom-left corner 

310 t 

311 | 

312 v 

313 a tmmA mm afl ■ alari Ire with data from the image 

314 ' IF mm fan* tmmU m* be loaded. THEN: 
313 / 

316: 4 L G+ h*mm mU** { A.ImAgcBitmap, pLine. iLineNum. IinageBiuiiap.B>4esPerLine)<0) 

317: ( 

318: // 

319: / rtim Error and abort 

320: V 

321: Error (MEM_ERK); 

322: L Fr»eBrtmap( AlmageBitmap y, 

323: drW plinc; 

324: THROW(theExoeption); 

323: // 

326: // END IF — late could not be loaded 

327: // 

328: } 

329: retura< pUne % 

330-7/ 

33 1 jl END - get a Imc from the Image 

332:// 
333:} 

334;// •••••••••••••••••••••••••• 

333 Jl Put a line m the Image 
3367/ •••••••••••••••••••••••••• 

337rvoid Image ::PutLtoe< BYTE •plragLine, UINT iLineNum ) 
338:{ 

339: CFikExooption *theExocotioo - new CFitcExceptioo; 

340: aaacrt( iUneNum < (UINT) lntageBhmap.Hcight \ 

341: // 

342: // IF the origin ia in the bottom left corner. THEN: 

343: // 

344: if( ImageBhraap. ViewPenpective «= BOTTOM LEFT ) 

343: ( 

346: // Adjust the row number based on the origin of the image 

347: iLineNum - ImagcBiunap.Height - iLineNum - 1 ; 

348: // 

349: //END IF — origin 01 bottom-left comer 

350: // 

351: } 

332: // 

353: // Put toe line into the imige 

334: //IF the line could not be placed, THEN: 

355: // 
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356: 
357 
358 
359 
360 
361 
362 
363 
364 
365 
366 
367 
368 
369:// 



370:// END - put a line in the Image 

371:// 

372:) 

373:// •••••••• 

3747/ Get a column from the Image 

3757/ 
376 



BYTE •lm*gc::GctCoiumn( UINT iColNum ) 



377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388: 

389 

390. 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410: 

411: 

412 

413 

414: 

415 

416 

417 

418 

419 

420: 

421 

422 

423 

424 

425 

426 



ifj L PulBhimp Row( AlmageBilmap. plmgUne. iUneNum, lma« C Bitmap. Bytes Per Une)<0) 
( 

// 

// Print Emir and abort 

// 

Error (M£M_ERR); 
L_FweBitmap( AUnageBitmap y, 
delete pLine; 

THROW( theExoeptioa y. 

II 

It END IF - line could not be placed 

// 
} 



UINT iRowNum; 
UINT tRow, 
BYTE pColPixel[3J; 
BYTE •pCon'mp; 

uiot( tColNum < (UINT) ImageBitmap. Width ); 

// 

// IF a column has not been allocated. THEN: 

// 

iT[ FlagAilocCol — FALSE ) 
{ 

// 

// Allocate the Column 

// 

prints •Allocating a column spaccta"); 
pCol - new BYTE(UnageBrtmap.Hei0it*31; 
FUgAllocCol - TRUE; 

// 

// END IF - a column has been allocated 

// 
) 

pCotTmp - pCol; 

// 

//LOOP through the number of pixeU in a apecujc column 

for( iRowNura-0; iRowNum < (UINT) InugeBrtmap.Hcight; iRowNunrM- ) 

iRow - iRowNum, 

// 

// IF the origin is m the bottom left comer, THEN: 

// 

ift ImageBmiup VicwPerspective — BOTTOM_LEFT ) 

II Adjust the row number baaed on the origin of the image 
iRow - ImageBhmap.Heigbt - 1 - iRowNum ; 

// 

// END IF - origin in bottom-left corner 

// 
) 

// 

// Read RGB pixel from the apecific cohimn 

// 

L_CktBitmapRowCoK Autu^cBitmap. pColPixel, iRow. iColNum, 3 ); 

// 

// Copy the .tGB pixel to the output buffer 
assert( pCoITmjvpCol < ImageBrtmap.Heighf3 ); 
•(pCoTTnm-M-) - •pColPixel 
•(pCoTTrop^) - •(pColPixel+1); 
•(pCorTmp++) - •(pColPixcl+2); 
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427: // END LOOP - number of pixel* in * specific column 
428: // 
429: ) 

430: n«unj(pCo4X 
43 VJl 

432:// END - get * column from (be Image 

433:// 

434:} 

435:// 

436V/ 

437V/ Put 4 column m the Image 

438 jf 

439:// 

440:void ImAge::PutColumn( BYTE* pCol UINT iCoINum ) 
441:{ 



442: UINT iRowNum; 

443: UINT iRow, 

444: BYTE -pColTmp; 

445: BYTE •pUncTmp; 

. aaaert( iCoINum < (UINT) ImageBhnup. Width ); 

447: assert* pCol !- NULL >. 

448: // 

449: // Create a tempormry line holder 

450: // 

45 1 : pUneTmp - no* BYTE(lmageBhrnap.BytesPerUne]; 

452: // 

453: // Copy the original column pointer 



454: // 

455: pCotTmp - pCol, 

456: // 

457: // LOOP through the number of pixels in a specific column 



458: // 

459: for{ iRowNum-O; iRowNum < (UINT) Image Bitmap. Height; iRowNum++ ) 
460: ( 

461: iRow - iRowNum; 

462: // 

463: //IF the origin is m the bottom left comer, THEN: 

464: // 

465: if( ImagcBctmap. ViewPcrtpeolivc — BOTTOM LEFT ) 

466: { 

467: // Adjust the row number based on the origin of the image 

468: iRow - ImageBifinap. Height * 1 - iRowNum ; 

469: // 

470: // END IF — origin in bottom-left corner 

471: // 

472: } 

473: // 

474: // Copy the RGB pixel to the output buffer 

475: // 

476: assect( pCofTmp-pCol < Image Bitmap. Height *3 ); 

477: // 

478: // Write the RGB pixel to the specific column 

479: // 

480: L OotBrtmapRov*( &IinageBctmap, pUncTmp, iRow, IroageBitmap. Bytes PcrLine X 

481: pUncTmp(0^iColNum # 3J - •pCoITmp++; 

482: pUneTmp! 1 +iCoINum*3J - •pCoITmp-M-; 

483: pLmeTmp(2-HColNum*3] - •pCotTmp++; 

484: L PutBhm*pRow( «fc Image Brtmap, pUncTmp, iRow, ImageBrtmap. Bytes PcrLine ), 

485: //" 

486: // NOTE : I suspect a bug in this Leadtools routine 

487: // 

488: //L PutBmnapRowCoK AlmagcBhrnap, Col Pixel iRow, iCoINura, 3 ); 
489: // 



490: // END LOOP - number of pixels in a specific column 
491: // 
492: } 

493: delete pLtneTmp; 
494:// 

495 Jl END -put a column in the Image 

496:// 

497:) 
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27/ Ococsi* Logic Classes 1.0 

3:// ... 
47/This»fltofda*a<* cmplomccts tome basic logic functions 

5 7/ for uae m the Genesis Warping Engine toftwarc 

67/ 

77/ Program Name: LOOICCPP 

87/ 

97/ Author James D. E. Ood 

107/ 

117/ - Initial release of the software. If you have any proWcrns, plcaae contact roc. 

127/ 

19930-^Microduphc. 20 0 Tc^n Centre Blvd., Suhe 400 

Markham. Ontario 

167/ 1311805 
yfj! (905)470-2742 

187/ 

197/ All righU reserved. 

207/ 

217/ Modifications: 

227/ 

23 M Aug. 8, 1994 • James Goel 

247/ -Initial release 

257/ 
26: 

27:tfndu<fc "logic.h" 
28:#includc <assert.h> 
29:#tnclude <stdio.h> 
30:#tnclude <afich> 
317/ 

327/ Create tome uacful type definition! 
337/ 

34 ^ypodcf unsigned char BYTE; 
33:typedef unsigned char PIXEL; 
36.1ypedef unsigned long ULONG; 
37*ypcdef unsigned a* WNT; 

• • — * 

427/ 

437/ Thai daa« emulate* a logic*! DFLOP. 

447/ 
45 . Jf — 

467/ D-Flop Clock transfers inputs to outputi 
477/ 

4«rvoid REG::Clk(wid) 
49:( 

50: Out -In; 
517/ 

527/ END - Constructor 

33:) 
54: 

557/ •••••••••••••••••••••••••••••••••••••••••••••• 

567/ Fiher Lookup Table 

577/ 
587/ 

597/ This cU^ iimui*tea the Filter Lookup tabic tn the warping 
607/ pateri. (reproduced below for ooropkscod) 

627/Tarlnc <Kt<-l I<»<-2 2<l<-4 4<t<-8 8<t<-l6 16<t<-32 

637/FUtSize 3 Tap 5 Tap 9 Tap 17 Tap 33 Tap 65 Tap 

647/FihFactor 0 12 3 4 5 

657/ 

667/ 1. Tarlnc it the apace between output Target pixels. 
677/ 2. FihSize if the FIR fiher tap size 

687/ 3. FihFactor is the number of right shifts required to divide the Tarlnc 

697/ inlo an appropriate Up Sample Increment. 

707/ 

7 1 7/ Constructor and Destructor 
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n-Jt 

73:LOOKUP::LOOKUP( void ) 

74:{ 

75:); 

76:LOOKUP::-UX)KUP( void ) 

77:{ 

78:}; 

79:// 

80:// Out — Output! a fitter factor based on the fTarlnc value 

81:// 

82:UINT LOOKUP::Out( void ) 
83:< 



84: aasert( fTarlnc > 0 X 

85: // 

86: // IF Tarlnc u in the fint category, THEN: 
87: // 

88. if{ (fTarlnc > 0) AA (fTarlnc <- 1)) 

89: // 

90: // Return t Fitter Factor of 0 for the 3 Tip FIR fitter 

91: // 

92: return 0; 

93: // 

94: // ELSE IF Tarlnc is in the second category, THEN: 

95: // 

96: else if ( (fTarlnc > 1) AA (fTarlnc <* 2)) 

97: // 

98: // Return a Fitter Factor of I for the 5 Tap FIR filter 

99: // 

100: return I; 

101: // 

102: // ELSE IF Tarlnc ts in the third category. THEN: 

103: // 

104: else if ( (fTarlnc > 2) AA (fTarlnc <- 4)) 

105: // 

106: // Return * Fitter Factor of 2 for the 9 Tap FIR fitter 

107: // 

108: return 2; 

109: // 

110: // ELSE IF Tarlnc is in the fourth category, THEN: 

111: // 

1 12: else if ( (fTarlnc > 4) AA (fTarlnc <- 8)) 
113: // 

114: // Return s Fitter Factor of 3 for the 17 Tap FIR filter 

115: // 

116: retum3; 

117: // 

118: //ELSE BF Tarlnc ism the fifth category, THEN: 

119: // 

120: else if ( (fTarlnc > 8) AA (fTarlnc <- 16)) 
121: // 

122: //Returns Fitter Factor of 4 for the 33 Tap FIR fitter 

123: // 

124: return 4; 

125: // 

126: // ELSE IF Tarlnc is in the sixth category, THEN: 

127: // 

128: else if ( (fTarlnc > 16) AA (fTarlnc <- 32)) 
129: // 

130: // Return s Filter Factor of 5 for the 65 Tap FIR filter 

131: // 

132: return 5. 

133: // 

134: // ELSE Tarlnc does not fit into any category, THEN: 

135: // 

136: else 

137: { 

138: // 

139: // Print internal error message 

140: // 

141: Jprintq stderr, "Fatal internal crror.Nn'); 

142: printfl; fTarlnc.. C/O]^, (Tarlnc y t 
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143: return 5; 

144: // 

145: // END IF - Filler Lookup Table 
146: // 
147: }; 
148:// 

149:// END FUNCTION- LOOKUP::Out 

150:// 

151:}; 

152: 

133:// 

154:// Tnree Input Comparator 
1557/ ••••••••••••••••••••••• 

1567/ 

1577/ Comparei three input* and output* the Urge* one. 
1587/ 

159:UINT TRICOMP::Out( void ) 
160: { 



161: UINT iTemp; 
162: // 

163: //Output the Urgest of three values 

164: // 

165: uTA>B) 

166: iTemp -A; 

167: eUe 

168: iTerop - B; 

169: 

170: if( iTemp > C) 

I7l : return iTemp; 

172: else 

173: return C; 



1747/ 

1757/ END FUNCTION - TRICOMP 

1767/ 

177:); 

178: 

1797/ •••*••••••**•••••*••*•• 

1807/ Two Input Comparator 
1817/ ••••••••••••••••••••••• 

1827/ 

1837/ Compare* two inpuU and output* the largest one. 
1847/ 

185:UINT BICOMP:Out( void ) 
186:( 

187: // 

188: //Output the largest of two values 

189: // 

190: if(A>B) 

191: return A; 

192: else 

193; return B; 

1947/ 

1957/ END FUNCTION - BICOMP 

1967/ 

197:); 

198: 

1997/ 

2007/ Upsampk Increment Generator 

2017/ 

2027/ 

203 7/ Right shifts Tarlnc by the Fitter Factor supplied by the 

2047/ Finer Lookup table. 

2057/ 

206:ULONO IncCen::Out( void ) 
207:{ 

208: ULONG lUpInc. 
209: // 

210: // Convert fTarlnc to an integer value for 
211: // faster snd essicr arthimetic caJoilsuocu 

lUplnc - (ULONG) ( fTarlnc • INDEX J>R£CISION_VALUE); 
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214: return ( IUptnc » FUtFsct ); 
2157/ 

2167/ END FUNCTION - Out 

2177/ 

218:}; 

219: 

220-7/ 

livj/hutpotau* 

2227/ ••••••••••••• 

2237/ 

2247/ Generates Interpolated pixels based oo the size of (he Fitter Factor 
225 Jl Generates a flag that indicates it baa oome to the end of a Taiinc 
2267/ value. 

, 2277/ 

228-7/ Default Constructor end Destructor 

2297/ 

230:lfdeipc4aie::IfiterpoUte< void ) 

231:{ 

232: // 

233: // Initialize private variables 

234: // 

233: Accum-OL; 
236: UpIncCount-O; 
237: NextPixPca-OU 
238: CurrPixPos-OL; 

2397/ 

2407/ END - Interpolate: :0)rtstructor 

2AVJI 

242:); 

243:Interpc4ate::MnterpoUte< void ) 

244:{ 

245:}; 

2467/ 

247:// Reset - Reset and initialize the intcrpolator 

2487/ 

249rvotd Interpolatc::Re*e^vo«d) 

230:{ 

231: // 

232: // Initialize private variables 

233: // 

254: Accum-OU 
235: UpIncCount-O; 
236: NcxtPixPos-OL; 
257: CurrPixPos-OL; 
238: iTarfncWidtn-O; 
239: // 

260: // Advance the current pixel value to be the next pixel value 
261: // 

262: SroCurrPix.Rod - *pSrcLaoe++; 
263: SrcCun-PucGm - •pSrcLine++; 
264: SrcCurrPix.Blu - «pSrcIjne++; 

265: // 

266: // Load SrcNextPix - GctNewSrc<pix) 
267: // 

268: SrcNextPix. Red - •pSrcUne++; 
269: SrcNextPix.Orn - •pSrclinc-J-f; 
270: SrcNextPix. Blu - •pSrcUne++; 
271: // 

272: // Initialize the Upaamplcd Pixel 

273: // 

274: tVHx.Red-SroCurrPix.Red; 
275: Upf4x.Gra-SrcCurrPix.Oni; 
276: UprHx.Bhi-SrcCurrPix.Blu; 

277:// 

2787/ END - InterpoUte:: Reset 

2797/ 

280:}; 

2817/ 

2827/ SetSccLine - Function to act the Source Line 

2837/ 

284 .-void ImerpoUle:^etSrc( BYTE •pLine ) 
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285:{ 
286: // 

287: // U>*d the source line into a private pointer 
288: // 

289: pSrcLine - pUnr, 
290: pSrcUncTmp - pUne; 
2917/ 

292*7/ END - u*crpouue::SctSrcUne 
2937/ 

294:); 

2957/ " ~ 

2967/CeoterTap - Indicate, when en interpolated pixel u the center Up for an 

2977/ FIR filter 

2987/ ■ ~ 

299:BOOLEAN u*errjoUle::CertkrTap( void ) 

300: ( 

302: // IF the interpolator tc creating the 6nt upearopled pixel OR 

303: // IF the interpolator is creating the last upaamplcd pixel THEN 

304: // 

305: if( UplncCount # 2 — (UINT) ( 2 « FihFact) ) 

306: { 

^g! // Initialize the UpSample mcremcnt counter and return * TRUE 

309! // 

3 10: UpIncCount - 0; 

311: return TRUE; 

312: ) 

313: eUe 

314: ( 

315: // 

316 //Advance the UptocCount 

317: If 

318: UpkcCount++; 

319: return FALSE; 

320: } 

321: // 

322: // END IF - interpolator 

323: // 

3247/ 

3257/ END - lnterpoUtor:CcnterTap 

3267/ 

327:); 

328: 

3297/ 

3307/ Out - Generates and up-aampk pixel 

3317/- 

332:Pixd lnterpoUte::Oonerale< void ) 

333:{ 

334: ULONQTeat; 

335: // 

336: // Accumulate the Up Sampled Incremenl 

337: // 

338: Accum 4- Uplnc; 

339' // 

340! //Theoext|Hxdpoeru«i«aoe^^ 

341: // 

342: NextPixPos - Gctlnt< Accum X 

344: // IF Next pixel position requires a new source pixel. AND 

345: //ffthencxtpixdhoi^a^paattheaidc/tte 

346* // 

347! tfl (NextPixPos > CurrPixPos) &A (NextPixPos < iTarlncWidth) ) 

348: { 

349: // 

35O: // Advance the current pixel position 

351: // 

332: CurrPixPos - NcxtPixPor, 

2U- //Advance the current pixel value to be the next pixd value 

355! // 
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356: SrcCurrPix.Red - SccNextPixRcd; 

357: SceCurrPix.Qrn - SrcNexlPixGra; 

358: SrcCurrPixBiu - SrcNextPixBhi; 

359: // 

360: II Load SrcNextPix - GctNcwSro(pix) 

361: // 

362: SrcNcxliHxRed - •pSrcLioc++; 

363: SrcNextPixGin - *pSecLtDe++; 

364: SrcNextPix Blu - *pSrcLinc^; 

365: // 

366: // END IF - NcxlPixPc* > CurrPixPo* 
367: // 
368: }; 
369: // 

370: // Calculate Interpolated pixel tutng linear mterpoUlion. 
371: //UprHx-SrcCurrPix* Fraction 
372: // 

373: Ten - GetFrac(Accum); 

374: UpPixRed - (BYTE) SrcCurrPixRed + (BYTE) 
'■■<>■•<* 375: (GetFi^A«^HSrcNexlPix.Red-SK PRECISION); 

376: UpPix Grn - (BYTE) SrcCurrPixGm ♦ (BYTE) 

377: (GetFrac(Accum)'(SrcNextPix.Gm - SrcCurrPixLGrn)>>INDEX PRECISION); 

378: UpPixBIu - (BYTE) SrcCurrPixBIu ♦ (BYTE) 

379: (Geu%ac<Accuni)»(SrcNexlPixBlu - SrcCurrPix.BIu>»INDEX^PRECISION); 

380: return UpPix; 

3817/ 

3827/ END FUNCTION - Imerpolalar.Out 

3837/ 

384:}; 

3857/ *********************** 
3867/ Left FIR Filter Kernel 
3877/ •••••••••••••**•••••••• 

3887/ 

3897/ Thk daa* cralea the left half of an Output Pixel baaed 
3907/ on the FIR filter indicated by the Filter Factor. 
3917/ 

3927/ Default Constructor and Deatnictor 

3937/ 

394:IxftJCcnieI::LeftKerncK void ) 

395: { 

396:}; 

397:LeftKeroel:-IuAKerncJ(void) 

398:( 

399:}; 

4007/ 

4017/ React - Initialize* the Left Kernel MAC 

4027/ 

" 403:void LeftKcrod::Raact< void ) 

-?7S 404: { 

405: // 

406: // Initialize private variable* 
407: // 

408: Double- FALSE; 
409: KcrnelPoe-0; 
410: AccumRedH); 
411: AooumOm-0; 
412: AccumBJu-0; 
4137/ 

4147/ END - LcftKerncl::React 

4157/ 

416:}; 

4177/ 

4187/ MAC - Filter a pixel 

4197/ 

420:void LeftKernel::MAC( Pixel UpPix ) 
42l:< 

422: UINTiEdgeFactor-1; 
423: // 

424: // IF the keracl u doubled. THEN: 
425: // 

426: ift Double— 'TRUE ) 
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427: // 

428: //"Double the UpSanipled Pixel 

429: H 

430: iEdgcF*c*or«2; 

431: // 

432: // END IF - double the kernel 

^33: // „ , 

4>4: // A ccumulate the FIR Coefficient • UpSanipled Pixel 

43tV /LouaJLed iLeftCocfl^hF*ctl(KenielP«] • UpPixRcd • iEdgeFatfor. 

437- A^umGru — fceftCod^FihF»ct](KewlPo«] • UpPix.Gra • iEdgeFactor. 

4 3g ; Aflc-nBhi — iUftCoefilFihFmrtJtKenielPocl • UpPix.BIu • iEdgeFactor, 

439: // 

440: //AJvwucio Ihe next kernel poahion 

441: // 

442: Mrs*<Pw~. 

444 // IF +*kMmi*mi*H* over. THEN: 



®0 445: / 

^ 446: ^ »Lm4r~ > (VTST) ( I « FUtFact)) 

A.,..:.; 44 7: 

44g. kir ai*at Israel 

449: 

450: kmm i*m - Q . 
451: /' 

452: // tKD IT -i«Hfc-f^4«va 

433: /' 
454:// 

455.// END - L«*k«t UAC 

456:// 
457:}; 
458://- 



459:// FiherPixH - ***** a FA* furi 

460://- 

46i:PixeI LenX*rw* FdurrW*>«d ) 

462:{ 

463- FihPix.lUd - (B^TE) ( AooumRed » UBYTE_MAX); 

464: FihPixBiu - (BYTE) < AccumBlu » UBYTE_MAX>, 

465: FihPixOni - (BYTE) (AccuroGm » UBYTE_MAX); 

466: // 

467: //CUmplheF»lM^p«uJtf required 

468: // 

469: aT AooumRad<0) 

470: FihPix. Rod - (BYTE) 0; 

471: if( AccuraOra<0) 

472: FthPuLGro - (BYTE) 0; 

473: AccumBlu <0) 

474: FihPixBlu - (BYTE) 0; 

475: 

476: ii(AocumRe<l> 65280) 

477: FihPix.Red - (BYTE) 255; 

478: tf( AccumGni > 65280 ) 

479: FUtPixGrn- (BYTE) 255; 

480: if( AcoumBlu > 65280) 

481: FihPix.Blu • (BYTE) 255; 

482: 

483: returo( FiKPix >. 
484:// 

485:// END - LeAKemel::FUtecPixel 

486:// 

487:}; 

488:// 

499 J/ Right FIR Filter Kernel 
490:// 

49VM , t 

492://Tnu das creaks the Right half of an Output Pixel owed 
493:// on the FIR 6her indicated by the Filter Factor. 
494:// 

495:// Default Cowtructor and Destructor 

496:// % 
4^:WfihtKeniel::RighlKenKl(^) SUBSTITUTE SHEET (RULE 26) 
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498: ( 
499:}; 

mftg!*KereI::~RighiKai>el( void ) 

50l:{ 

302:}; 

503:// 

504:// React - Initialize the Right Kernel MAC 

505:// 

506.W1 RightKeroel::ReMt< void ) 
507:{ 

508: // 

509: // Initialize private variable* 
510: // 

511: Double -FALSE; 
512: KemelPor-0; 
513: AooumRod-0; 
514: AooumOra-0; 
515: AooumBIu-O. 
516:// 

5177/ END - RighUCeroel::Reset 

5187/ 

319:}; 

3207/ 

521:// MAC - Filler a pixel 

322:// 

523:void RightKemel::MAC( Pixel UpPix ) 
524:( 



525: UINT iEdgcFactor - 1; 

326: // 

527: //IF the kernel is doubled, THEN: 

328: // 

529: if( Double— TRUE ) 
330: // 

531: //Double the UpSampled Pixel 

532: // 

533: iEdgeFactor-2; 

534: // 

535: // END IF — double the kernel 



336: // 

537: // Accumulate the FIR Coefficient * UpSampled Pixel 

538: // 

539: AceumRcd +- u^ghtCocflIFihTact][KcriieiPo.] • UpPixRed • iEdgeFacton 

540: AccumGrn -h- iRightCocff|Filtf *ct][Korr>ciPo«] • UpPix-Gm • iEdgcFactor. 

541: AooumBlu +- iRightCocfl[FihFactJ[KernelPo.] • UpPix_Blu • iEdgcFactor; 

542: KernelPorH-; 

343: // 

544: // IF the kernel baa rolled over. THEN: 

545: // 

346: il(KecnelPo«>(UINT)(l<<FihFact)) 



547: // 
348: //Initialize the Kernel 

349: // 
550: KernelPor-O; 

551: // 

532: //END IF -kernel has rolled over 

553: // 



354:// 

555:// END - RighlKernel::MAC 

556:// 

337:}; 

358:// 

5597/ FtherPixel - Return a Finer Pixel 

5607/ 

561:Pixcl Ri«hlKernel::FUterPixel( void ) 
562: ( 



563: // 

564: // Convert the accumulated filtered pixel to a byte 

565: // 

566: FihPix.Rcd - (BYTE) (AceumRcd » UBYTE_MAX); 

567; FiKPix.Grn - (BYTE) (AccumOrn » UBYTE MAX); 

568: FihPix.Blu - (BYTE) (AceumBhi » UBYTE MAX); 
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569: // . J 

570: // CUn*> the FiHorod pixsl if rotfiircd 

571: // 

572: if( AfloumRed < 0) 

373: FiltPix.Rod - (BYTE) 0; 

574: AocumOni < 0) 

575: FihPixOni - (BYTE) 0; 

576: LfC A^cumBha < 0) 

577; FihPixBlu - (BYTE) 0; 

578: 

579: if(AoeumRe<l> 65280) 

580: FUtPix.Red • (BYTE) 255; 

581; ifl[AocumQtTi> 65280) 

582: FiltPixGni- (BYTE) 255; 

583: if( AceumBlu > 65280) 

584; FtltPixBlu - (BYTE) 255; 

585: retting FUlPix); 

586:// 

587:// END - Rigt*Kcmcl::FiHcrPixc! 

5887/ 

589:); 
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I'M 

VJI Genesis Logic CUmc* 1.0 

yji 

A M This set of cUmct implement! tome basic logic functions 
$ jl for use in the Genesis Warping Engine software. 

6V/ 

7:// Program Name: LOOICCPP 
8:// 

9:// Author Jubm D. E. God 
107/ 

ll'M - Initial release of the software. If you have any problems. ple*ie contact mc. 

u-y/ 

Uy/CopyKfnt(c) 1993 Genesis Microchip Inc. 

200 Town Centre BIvdL, Suite 400 
IW Markhara, Ontario 

L3R80S 
(905)470-2742 

It/' Al nfm tmwnm& 

»Jt 

aiy/Hiif ii m 

22-tf 

2Jr t. |*M - James Gocl 

24 y/ - Initial release 

23 Jt 

24 

27 *aWrf_LDGICCUaS H_ 
21 MaW _LOGICCl-AK H_ 
29< 

30: 
31 1 

32 *«m«i Movarf fr —lunn Clasaes require C++ compilation (use a .cpp suffix) 
34: 

33aypa«Wiiii piiti r BYTE; 
3<S.Typ«drf— iiariM ofaar PIXEL; 
37.typ.aW. iisji I long ULONG; 
38typ«oW ,„■«»■ 1 «m UINT; 
39: 
40// 

4 1 JJ Define a Piad «wait 

42:// 

43 typodef struct puoal tag 

44:( 

45: BYTE Red; 
46: BYTE Gro. 
47: BYTE Blu; 
48:} Pixel; 
49: 

5o-y/ 

31// Setup a boolean type with TRUE-1 and FALSE-0 
32// 

33.typodef unsigned mt BOOLEAN; 
54:#definc TRUE 1 
33:#define FALSE 0 
36: 

37y/ 

58// Define the fractional prccuioo of Tsrlnc 

39// 

60:#o^fineINDEXJ>R£CISION 16L // 16 fractional bhi of precision 

6l:#define INDEX_PR£CISION VALUE (1L « INDEX PRECISION) 

62:tfdefinc INDEX PRECISION" MASK ((1L « INDEX PRECISION)- I L) 
63: 

64y/ 

63 // Define the number of bits per byte 
66// 

67:#dcfine UBYTE MAX 8 
68: 

69:«define GoUnt(Vaiue) ((ULONGX(Value) » INDEX_PRECISION)) 

70:#oVfineOetfrac(Value) ((Value) A INDEX PRECISION MASK) 

71: ~ 
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72:// 

73:// xl FIR Fiher Cocfficionts 
74:// 

75:r¥define MAXFILTER 6 // Maximum filter type number 

76 #define MAX TAPS 33 // Maximum number of filter Up* (KenuU doubled) 

77: 
787/ 

797/ Left Kernel Coefficient* 

80:// 

8 1 :const long U>ftCoefi[MAX_FILTER][MAX_TAPSJ - 
82:{ 

83: {0. 128, 0. 0, 0. 0, 0, 0, 0. 0, 0. 0, 0. 0, 0, 0, 0, 0. 0, 0, 0. 0, 0, 0. 0, 0, 0, 0, 0, 0. 0, 0. 0}, 
84: {-16. 64*, 80, 0. 0, 0, 0, 0. 0, 0. 0, 0. 0, 0, 0, 0, 0, 0, 0, 0. 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), 
85: { -7. 3. 30, 63, 39. 0. 0. 0. 0, 0. 0. 0, 0, 0, 0, 0, 0, 0. 0, 0, 0, 0, 0, 0, 0. 0, 0. 0, 0. 0. 0, 0, 0), 
86: { -3. -2, U 7, 15. 23, 3 1. 37. 19. 0. 0. 0. 0. 0. 0, 0, 0. 0. 0. 0, 0. 0. 0, 0. 0, 0. 0. 0, 0. 0, 0. 0. 0}. 
87 {-1,-1,-1. 0. 0. 2. 3. 5. 7,9,11,13,15,163344.9.0,0,0,0, 0. 0.0.0.0.0,0,0,0.0,0,0), 
88: { 1. 1. 1. 2. 2, 3. 3. 4, 4, 5. 5, 6, 6. 7. 7. 8, 8. 8. 8. 9. 9. 16. 5. 0. 0. 0. 0, 0, 0. 0. 0. 0. 0} 
89:}; 
907/ 

91:// Right ICemcl Coefficient* 

92:// 

93:eomt long i Ri ghtCoeflJMAX_FILTE R] [ MAX T APS ] - 
94:{ 

95: {128. 0. 0, 0, 0, 0, 0, 0, 0, 0. 0, 0, 0, 0. 0, 0, 0, 0, 0, 0. 0. 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0. 0>, 
96: {80. 64^ -16, 0, 0, 0. 0, 0, 0, 0, 0, 0. 0, 0. 0, 0, 0. 0, 0, 0, 0, 0, 0, 0. 0, 0, 0, 0, 0, 0. 0. 0. 0), 
97: {39^ 63. 30. 3. -7, 0, 0, 0, 0, 0. 0. 0. 0, 0, 0, 0, 0, 0, 0. 0. 0, 0. 0, 0. 0. 0, 0. 0. 0. 0, 0. 0. 0), 
98: {19^ 37. 31.23. 15, 7. 1,-2,-3,0,0,0.0.0,0,0,0.0.0.0.0,0.0,0,0.0.0.0.0.0.0.0.0}. 
99: {9. 24. 33, 16. 15. 13. U. 9. 7. 5, 3, 2. 0. 0.-1.-1.-1. 0, 0, 0. 0, 0. 0, 0. 0, 0. 0. 0. 0. 0. 0, 0, 0}, 
100: {5. 16. 9, 9. 8. 8. 8. 8. 7. 7. 6. 6, 5, 5, 4, 4, 3, 3, 2. 2, 1, 1, 1, 0. 0, 0. 0, 0. 0. 0. 0. 0, 0} 
101:}; 
102: 

103V/***— ••••—*••••••••••••••••••••••*•••••••*••••••••••••••••••*••••••••*•• 

104:// D-Flop Register Class 

106:// 

1077/ This das* simulate* a logical DFLOP. 
108:// 

109:cU*aR£G 
110:{ 

lllrprivate: 
112:puWic: 
113: // 

114: // Input and Output Register Value 
115: // 

116: UINT In; 
117: UINT Out; 
118: // 

119: // Class Constructor and Destructor 
120: // 

121: REG(voidX 
122: -R£0( void); 
123: // 

124: // Clock register value 
125: // 

126: void Clk( void); 
127*// 

1287/ END CLASS - REG 

1297/ 

130:}; 

131: 

132:// *************************************************************************** 

1337/ Filter Lookup Table 

135-// 

1367/ This class simulates the Filter Lookup table in the warping 

1377/ patent (reproduced below for completeness) 

1387/ 

1397/TarInc 0«=l I<t<«"2 2<l<«=4 4<t<«8 8<t<-16 16<l<-32 
14077 FiitSize 3 Tap 5 Tap 9 Tap 17 Tap 33 Tap 65 Tap 
1417/FihFactor 0 1 2 3 4 5 

1427/ 
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143:// 1. Tarlnc b the spaoe between output Target pixels. 
1447/2. FiltSixe b the FIR fiher Up size 

145-7/ 3. FtftF actor is the number of right shifts required to divide the Tarlnc 

1467/ into an appropriate Up Sample focranent 

1477/ 

148:da« LOOKUP 
149:{ 

150:privmle: 
15L-public: 



152: 


// 


153: 


// Cocvtructor and Destructor 


154: 


// 


155: 


LOOKUP(void); 


156: 


~LOOKUP(void> f 


157: 


// 


158: 


// Output Target pixel spacing 


159: 


// 


160: 


float fTarlnc; 


161: 


// 


162: 


// Mapped Fitter Factor for the input Tarlnc 


163: 


// 


164: 


UINT Out(void); 


165:); 




166: 





1677/ 

1687/ Three Input Comparator 
1697/ ••••••••••••••••••••••• 

1707/ 

1717/ Compare* throe inputs and output* the Largest one. 
1727/ 

173xlaa* TRJCOMP 
174:{ 

175 private: 
176:puWic: 
177: // 

178: // Constructor and Destructor 
179: // 

180: TRJCOMr\void); 
181: -TRICOMP(void); 
182: // 

183: //Three input variables 
184: // 

185: UDMTAAC. 
U v';j 186*: // 

187: //Generate the output 
118: // 

189: UINT Out( void); 
M 190:}; 

£33 191: 

1927/ 

1937/ Two Input Comparator 

1947/ 

1957/ 

1967/ Compares two inputs and outputs the largest one. 
1977/ 

198:cUss BICOMP 
199:{ 

200:orivate: 
201:pubuc: 



202: 


// 


203: 


//Constructor and Destructor 


204: 


// 


205: 


BICOMFXvoid); 


206: 


~BICOMP(void); 


207: 


// 


208: 


//Two input variables 


209: 


// 


210: 


UINT A3; 


211: 


// 


212: 


// Generate the output 


213: 


// 
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214: UINTOut(voi<l); 

215:}; 

216: 

217:// ••••••••« 



218://Upsamplcl 
219:// •••••••••••*•••••••••••••••*• 

220:// 

22 VJt Right shifts Tarfnc by the Filter Factor supptied by the 
2227/ Filter Lookup tsble. 

223:// 

224:cUm IncGen 
225:{ 

226:private: 
227:pubtic: 
228: // 

229: //t^saiuple Filter Factor 
230: // 

231: UINT FUtFact; 
232: float flarlnc, 
233: // 

234: // Output in Upumpte Increment 
235: // 

236 : ULONG Out( void ); 

237:); 

238:// 

239:// Interpolator 
240^/ ••••*••••••*• 

241:// 

242 ji Generates Interpolated pixels based on the tize of the Filter Factor 
243:// Generates a flag that indicates it has come to the end of a Tarlnc 
244:// value. 
245:// 

246xUss Interpolate 
247:{ 

248:private: 

249: UINT UplncCount; 

250: ULONG CurrPixPoa, NextPixPor, 

251: Pixel SrcCuirPix.SrcNexlPix; 

252: BYTE •pSrcUne. •pSrclineTmp; 

253:puMk: 

254: // 

255: // Upsampie Fther Factor 

256: // 

257: ULONG Accum; 
258: UINT FihFact; 
259: float fTarlnc. 
260: Pixel UpPix; 
261: BYTE bPixd; 
262: ULONO Uplnc, 
263: UINT iTsrInc Width; 
264: // 

265: //Constructor and Destructor 

266: // 

267: Interpc4ate<void); 
268: -toterpoUle<voidX 
269: // 

270: //Rea^UtheinterpoUlorloaknoMmsiate 
271: // 

272: void He*et( void); 
273: // 

274: //FiinctiontoactthcpouilertotheSou 

275: // 

276: void SctSrc( BYTE •pUnc y. 
277: // 

278: // Generate an Interpolated Pixel 
279: // 

280: Pixel Gencrate< void ); 
281: // 

282: // Flag indicates when the center kernel 
283: // Up ixaerpoUled pixel is bong gcoersled 
284: // 
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285: BOOLEAN CeoterTap( void y t 
296:): 

287// *************************** 

288// Left FIR Filter Kernel MAC 
289// *************************** 

290:// 

29 1 // Thii claa* cneatee the left half of an Output Pixel baaed 
292// an the FIR filter indicated by the Filter Factor 
293// 

294:cUaa LeftKeroel 

295;{ 

296 .-private: 

297: UINT KerodPoc, 
298.-puWic: 

299: long AooumRed; 

300: long AecumOrn; 

301: long AccumBlu; 

302: BOOLEAN Double. 

303: Pixel FihPix; 

304: UINT FihFact; 

305: // 

306: // Cooatnictor and Dcatmctor 

307: // 

308: LeftKerneK void ); 
309: ~LeftKrroel< void); 
310: // 

311: // React MAC to a known atate 

312: // 

313: void Reaet( void); 
314: // 

315: // Generate a filtered pixel 
316: // 

317: void MAC( Pixel UpPix); 
318: // 

3 19: // Return a Filtered Pixel 
320: // 

321: Pixel Filter Pixel( void ); 
322// 

323// END CLASS - LeflKernclMAC 

324// 

325:); 

326// **********•**•***••**•*•••• 

327// Right FIR Filter Kernel MAC 
32gy/ •*••«•«•*••••••••*•*«•«•*•* 

329// 

330// Tbia ciaaa crealee the right half of an Output Pixel baaed 
331// on the FIR filter indicated by the Fiher Factor. 
332// 

333:cUaa RigbUCemel 
334:{ 

335.-prrvatc: 

336: UINT KeroclPor. 
337,-pubUc: 

338: long AooumRed; 

339: long AccumGra; 

340: long AcoumBIu; 

341: BOOLEAN Double; 

342: Pixel FthPix; 

343: UINT FUtFact, 

344: // 

345: // Constructor mud Deatructor 
346: // 

347: RightKert>eJ<void); 
348: -4UghtKerocl( void); 
349: // 

350: //React MAC to a known aute 

351: // 

352: void Reaet( void); 

353: // 

354: // Ocnorate a filtered pixel 

355: // 
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356: void MAC( Pixel UpPix); 

357: // 

358: // Return a Filtered Pixel 

359: // 

360: Pixel FiherPixel( void); 
361:// 

162'Jf END CLASS - RigfatKmeiMAC 

363:// 

364:>; 
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17/ 

2M ImJLgc CUs* - Encapsulates the Leadtools Image Functions 

37/ 

47/7hia image daas greatly improves the LeadTooIs user interface by 
37/ providing a simplified and consistent set of C++ member function* 
67/ for manipulating images, 

77/ 

87/ Program Name: IMAGE.CPP 

97/ 

^ ? 107/ Author: James D. E> Ood 

117/ 

127/ - Initial release of the software. If you have any problem*, please oonUct me. 

137/ 

1 47/ Copyright (c) 1993 Genesis Microchip Inc. 

137/ 200 Town Centre Blvd. Suite 400 

167/ Markham. Ontario 

177/ L3R8G5 

187/ (905)470-2742 

197/ 

207/ All rights reserved 

v- V-ii 217/ 

@i 227/M«itfc*4iom: 
^ 237/ 

247/ Aug. 5,1994 - James Goel 

257/ - initial release 

267/ 

277/ Ciaas definition 

287/ 
297/ 

30:*defme FOR DOS 
3t:#o*dude-1_tooUp.ir 
32:#mdu4kn^bftmapJi- 
33:#inciuoe "I errorii- 

347/ 

337/ Setup LEAD TOOLS for DOS definitions 

367/ 

37:#imdef_IMAQECL7\SS H__ 
38:#<fafjne __IMAGECLASS~H_ 
39;#endtf 
40: 

4I:#imdef_c*^uaplus 

42:#error Microsoft Foundation Classes require C++ compilation (use a .cop suffix) 

43:«en£f 

44: 

43*ypedef unsigned char BYTE; 
46*ypodef unsigned char PIXEL; 
47^ypedef unsigned long ULONO; 
48*ypedcf unsigned eat UINT; 
49*ypodef unsigned inl BOOLEAN; 
50: 

51:#defineTRUEl 
52:*define FALSE 0 



53: 






54:cU 


as Image 




53:( 






56:prtvate: 




37: 


BITMAPHANDLE IniageBrtmap; 


58: 


char 


txhInputFdename(l3J; 


59: 


char 


pchOutputFUeoame[l3); 


60: 


UINT 


iWidtfa; 


61: 


UINT 


iHcight; 


62: 


UINT 


iline; 


63: 


BYTE 


•pUne; 


64: 


BYTE 


•pCcd; 


65: 


BYTE 


•pTGALme; 


66: 


FILE 


VT^rgt; 


67: 


BOOLEAN FlagAllocLinc; 


68: 


BOOLEAN FlagAllocCol; 


69: 


void 


Error( mt ermum X 


70: 







71 public: 
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72: // 

73: //CUm Constructor and Deatructor 

74: // " 

75: lmage(void); 

7$: -*mage<void); 

77: // 

78: //Iw allocated Bitmap and memory 

79: // 

SO: void FreeMeraory(void); 

81: // 

82: //Load end Save image file 

83: // 

84: void Load( char *pcninput ); 

85: void Save< char V****** 

86: // 

87: //Get image Width and Height 

88: // 

89: UINT WidlKvotdX 

90: UINT Height void y, 

91: // 

92: // Oct and Put image lines 

93: // 

94: BYTE *GetUne( UINT UineNum ); 

95: void PutLinc( BYTE •pLine, UINT iLmeNum ); 

96: // 

97: //Get and Put image columns 

98: // 

99: BYTE •GetCc4umn( UINT iUneNum ); 

100: void PutColumn( BYTE •pUne, UINT iUncNura ); 

101: // 

102: //Create an empty image 

103: // 

104: void Oreaie( UINT iWidth, UINT iHcight X 

105-y/ 

106:// END CLASS -I 
107:// 
108:}; 
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1. A sample sequence warping system for receiving 

and warping an input sequence of source samples to 
produce an output sequence of target samples, comprising: 

a) first means for receiving successive target 
increments representing desired target spacing between 
respective ones of said target samples , and in response 
generating successive filter factor values representing 
respective FIR filters having sufficient filter 
coefficients to generate respective ones of said target 
samples with said desired target spacing; 

b) second means for receiving said successive 
filter factor values and in response selecting 
predetermined left and right half kernels of said 
respective FIR filters for generating said respective 
ones of said target samples; 

c) third means for receiving said successive 
target increments and said successive filter factor 
values and in response generating a predetermined number 
of upsampled intermediate samples from said input 
sequence of source samples; and 

d) fourth means for applying said upsampled 
intermediate samples to said predetermined left and right 
half kernels of said respective FIR filters for 
generating said output sequence of target samples. 

2. The sample sequence warping system of claim 1, 
wherein said first means comprises a filter look-up table 
for correlating each said target increment with a 
respective one of said filter factor values and sizes of 
said FIR filters. 

3. The sample sequence warping system of claim 2, 
wherein said filter look-up table generates a filter 
factor value of 0 when said target increment is between 0 
and 1. 
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4 # The sample sequence warping system of claim 2, 

wherein said filter look-up table generates a filter 
factor value of 1 when said target increment is between 1 
and 2. 

5 

5> The sample sequence warping system of claim 2, 

wherein said filter look-up table generates a filter 
factor value of 2 when said target increment is between 2 
and 4. 

10 

6. The sample sequence warping system of claim 2, 
wherein said filter look-up table generates a filter 
factor value of 3 when said target increment is between 4 
and 8. 

15 

7. The sample sequence warping system of claim 2, 
wherein said filter look-up table generates a filter 
factor value of 4 when said target increment is between 8 
and 16. 

20 

8. The sample sequence warping system of claim 2, 
wherein said filter look-up table generates a filter 
factor value of 5 when said target increment is between 
16 and 32. 

25 

9. The sample sequence warping system of claim 1, 
wherein said second means further comprises means for 
selecting said predetermined left and right half kernels 
of the largest one of said respective FIR filters on 

30 either side of each of said target samples. 

10. The sample sequence warping system of claim 1, 
wherein said second means further comprises means for 
selecting said predetermined left and right half kernels 

35 from respective ones of said FIR filters on either side 
of each of said target samples. 
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11. The sample sequence warping system of claim 9, 
wherein said means for selecting said predetermined left 
and right half kernels further comprises: 

i) first register means connected to said first 
means, for receiving and delaying said successive filter 
factor values by one sample period; and 

ii) first comparator means connected to said 
first means and said first register means, for selecting 
the largest of two successive filter factor values output 
from said first means and said first register means, and 
in response outputing successive largest ones of two said 
successive filter factor values. 

12. The sample sequence warping system of claim 11, 
wherein said third means further comprises: 

iii) second register means connected to said 
first register means for receiving and delaying said 
successive filter factor values by a further sample 
period; 

iv) second comparator means connected to said 
first means, said first register means and said second 
register means, for selecting the largest of three 
successive filter factor values output from respective 
ones of said first means, said first register means and 
said second register means, and in response outputing 
successive largest ones of three said successive filter 
factor values; 

v) an upsample increment generator for 
receiving said successive target increments and said 
successive largest ones of three said successive filter 
factor values and in response right shifting said 
successive target increments respectively by said 
successive largest ones of three said filter factor 
values to form successive upsample increment values; and 

vi) interpolator means for receiving said 
successive binary target increments, said successive 
upsample increment values and said input sequence of 
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source samples, and in response interpolating said input 
sequence of source samples to produce said predetermined 
number of upsampled intermediate samples spaced apart by 
a respective one of said upsample increment values, said 
5 predetermined number being equivalent to a respective one 
3 of said target increments divided by said respective one 

of said upsample increment values. 

13 m The sample sequence warping system of claim 12, 

10 wherein said fourth means further comprises: 

vii) third register means connected to said 
second comparator means for receiving and delaying said 
successive largest ones of two said successive filter 
factor values by one sample period; 
15 viii) fourth register means connected to said 

third register means for receiving and delaying said 
successive largest ones of two said successive filter 
factor values by a further sample period; 

ix) first half digital filter kernel multiplier 
20 accumulator means connected to said third register means 

and said interpolator means, for receiving and 
multiplying a first plurality of said upsampled 
intermediate samples by successive ones of said filter 
m: coefficients identified by said successive largest ones 

25 of two said successive filter factor values delayed by 

one sample period and in response generating a first 
plurality of intermediate product values , summing said 
first plurality of intermediate product values and in 
response generating a succession of accumulated 

30 intermediate product values; and 

x) second half digital filter kernel multiplier 
accumulator means connected to said first half digital 
filter kernel multiplier accumulator means, said fourth 
register means and said interpolator means, for receiving 

35 and multiplying a second plurality of said successive 

upsampled intermediate samples by successive ones of said 
filter coefficients identified by said successive largest 
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ones of two said successive filter factor values delayed 
by a further sample period and in response generating a 
second plurality of intermediate product values, summing 
said second plurality of intermediate product values and 
in response generating a further succession of 
accumulated intermediate product values, and summing said 
succession of accumulated intermediate product values and 
said further succession of accumulated intermediate 
product values to generate said output sequence of target 
samples. 

14. The sample sequence warping system of claim 10, 
wherein said third means further comprises: 

i) an upsample increment generator for receiving 
said successive target increments and said successive 
filter factor values and in response right shifting said 
successive target increments respectively by said 
successive filter factor values to form successive 
upsample increment values; and 

ii) interpolator means for receiving said successive 
target increments, said successive upsample increment 
values and said input sequence of source samples, and in 
response interpolating said input sequence of source 
samples to produce said predetermined number of upsampled 
intermediate samples spaced apart by a respective one of 
said upsample increment values, said predetermined number 
being equivalent to a respective one of said target 
increments divided by said respective one of said 
upsample increment values. 

15. The sample sequence warping system of claim 14, 
wherein said fourth means further comprises: 

iii) register means connected to said first 
means, for receiving and delaying said successive filter 
factor values by one sample period; 

iv) first half digital filter kernel multiplier 
accumulator means connected to said register means and 
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said interpolator means, for receiving and multiplying a 
first plurality of said upsampled intermediate samples by 
successive ones of said filter coefficients identified by 
said successive filter factor values delayed by one 
5 sample period and in response generating a first 

plurality of intermediate product values, summing said 
first plurality of intermediate product values and in 
response generating a succession of accumulated 
intermediate product values; and 

10 V ) second half digital filter kernel multiplier 

accumulator means connected to said first half digital 
filter kernel multiplier accumulator means, said register 
means and said interpolator means, for receiving and 
multiplying a second plurality of said successive 

15 upsampled intermediate samples by successive ones of said 

filter coefficients identified by said successive filter 
factor values delayed by said one sample period and in 
response generating a second plurality of intermediate 
product values, summing said second plurality of 

20 intermediate product values and in response generating a 

further succession of accumulated intermediate product 
values, and summing said succession of accumulated 
intermediate product values and said further succession 
of accumulated intermediate product values to generate 

25 said output sequence of target samples. 

16 . a sample sequence warping method for receiving 

and warping an input sequence of source samples to 
produce an output sequence of target samples, comprising 

30 the steps of: 

a) receiving successive target increments 
representing desired target spacing between respective 
ones of said target samples, and in response generating 
successive filter factor values representing respective 

35 FIR filters having sufficient filter coefficients to 

generate respective ones of said target samples with said 
desired target spacing; 
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b) receiving said successive filter factor 
values and in response selecting predetermined left and 
right half kernels of said respective FIR filters for 
generating said respective ones of said target samples; 

c) receiving said successive target increments 
and said successive filter factor values and in response 
generating a predetermined number of upsampled 
intermediate samples from said input sequence of source 
samples; and 

d) applying said upsampled intermediate samples 
to said predetermined left and right half kernels of said 
respective FIR filters for generating said output 
sequence of target samples. 

17. The sample sequence warping method of claim 16, 
wherein said step of receiving said successive filter 
factor values and in response selecting predetermined 
left and right half kernels of said respective FIR 
filters for generating said respective ones of said 
target samples further comprises selecting said 
predetermined left and right half kernels of the largest 
one of said respective FIR filters on either side of each 
of said target samples. 

18. The sample sequence warping method of claim 16, 
wherein said step of receiving said successive filter 
factor values and in response selecting predetermined 
left and right half kernels of said respective FIR 
filters for generating said respective ones of said 
target samples further comprises selecting said 
predetermined left and right half kernels from respective 
ones of said FIR filters on either side of each of said 
target samples. 

19. The sample sequence warping method of claim 17 , 
wherein said step of selecting said predetermined left 
and right half kernels further comprises: 
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i) receiving and delaying said successive 
filter factor values by one sample period; and 

ii) selecting the largest of two successive 
filter factor values and in response outputing successive 

5 largest ones of two said successive filter factor values. 

20. The sample sequence warping method of claim 19, 

wherein said step of receiving said successive target 
increments and said successive filter factor values and 
10 in response generating a predetermined number of 

upsampled intermediate samples from said input sequence 
of source samples further comprises: 

iii) receiving and delaying said successive 
filter factor values by a further sample period; 

15 iv) selecting the largest of three successive 

filter factor values and in response outputing successive 
largest ones of three said successive filter factor 
values; 

v) receiving said successive target increments 
20 and said successive largest ones of three said successive 

filter factor values and in response right shifting said 
successive target increments respectively by said 
successive largest ones of three said filter factor 
values to form successive upsample increment values; and 
25 vi) receiving said successive binary target 

increments, said successive upsample increment values and 
said input sequence of source samples, and in response 
interpolating said input sequence of source samples to 
produce said predetermined number of upsampled 
30 intermediate samples spaced apart by a respective one of 

said upsample increment values, said predetermined number 
being equivalent to a respective one of said target 
increments divided by said respective one of said 
upsample increment values. 
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21. The sample sequence warping method of claim 20, 

wherein said step of applying said upsampled intermediate 
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samples to said predetermined left and right half kernels 
of said respective FIR filters for generating said output 
sequence of target samples further comprises: 

vii) receiving and delaying said successive 
largest ones of two said successive filter factor values 
by one sample period; 

viii) receiving and delaying said successive 
largest ones of two said successive filter factor values 
by a further sample period; 

ix) receiving and multiplying a first plurality 
of said upsampled intermediate samples by successive ones 
of said filter coefficients identified by said successive 
largest ones of two said successive filter factor values 
delayed by one sample period and in response generating a 
first plurality of intermediate product values, summing 
said first plurality of intermediate product values and 
in response generating a succession of accumulated 
intermediate product values; and 

x) receiving and multiplying a second plurality 
of said successive upsampled intermediate samples by 
successive ones of said filter coefficients identified by 
said successive largest ones of two said successive 
filter factor values delayed by a further sample period 
and in response generating a second plurality of 
intermediate product values, summing said second 
plurality of intermediate product values and in response 
generating a further succession of accumulated 
intermediate product values, and summing said succession 
of accumulated intermediate product values and said 
further succession of accumulated intermediate product 
values to generate said output sequence of target 
samples . 

22 • The sample sequence warping method of claim 18, 

wherein said step of receiving said successive target 
increments and said successive filter factor values and 
in response generating a predetermined number of 
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upsampled intermediate samples from said input sequence 
of source samples further comprises: 

i) receiving said successive target increments 
and said successive filter factor values and in response 
right shifting said successive target increments 
respectively by said successive filter factor values to 
form successive upsample increment values; and 

ii) receiving said successive target 
increments, said successive upsample increment values and 
said input sequence of source samples, and in response 
interpolating said input sequence of source samples to 
produce said predetermined number of upsampled 
intermediate samples spaced apart by a respective one of 
said upsample increment values, said predetermined number 
being equivalent to a respective one of said target 
increments divided by said respective one of said 
upsample increment values. 

23. The sample sequence warping method of claim 22, 

wherein said step of applying said upsampled intermediate 
samples to said predetermined left and right half kernels 
of said respective FIR filters for generating said output 
sequence of target samples further comprises: 

iii) receiving and delaying said successive 
filter factor values by one sample period; 

iv) receiving and multiplying a first plurality 
of said upsampled intermediate samples by successive ones 
of said filter coefficients identified by said successive 
filter factor values delayed by one sample period and in 
response generating a first plurality of intermediate 
product values, summing said first plurality of 
intermediate product values and in response generating a 
succession of accumulated intermediate product values; 
and 

v) receiving and multiplying a second plurality 
of said successive upsampled intermediate samples by 
successive ones of said filter coefficients identified by 
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said successive filter factor values delayed by said one 
sample period and in response generating a second 
plurality of intermediate product values, summing said 
second plurality of intermediate product values and in 
response generating a further succession of accumulated 
intermediate product values, and summing said succession 
of accumulated intermediate product values and said 
further succession of accumulated intermediate product 
values to generate said output sequence of target 
samples. 
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